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Fractional White-Noise Limit and Paraxial Approximation for 

Waves in Random Media 

Christophe Gomez* and Olivier Pinaud^ 


Abstract 

This work is devoted to the asymptotic analysis of high frequency wave propagation 
in random media with long-range dependence. We are interested in two asymptotic 
regimes, that we investigate simultaneously: the paraxial approximation, where the wave 
is collimated and propagates along a privileged direction of propagation, and the white- 
noise limit, where random fluctuations in the background are well approximated in a 
statistical sense by a fractional white noise. The fractional nature of the fluctuations is 
reminiscent of the long-range correlations in the underlying random medium. A typical 
physical setting is laser beam propagation in turbulent atmosphere. Starting from the 
high frequency wave equation with fast non-Gaussian random oscillations in the velocity 
field, we derive the fractional Ito-Schrodinger equation, that is a Schrodinger equation 
with potential equal to a fractional white noise. The proof involves a fine analysis of 
the backscattering and of the coupling between the propagating and evanescent modes. 
Because of the long-range dependence, classical diffusion-approximation theorems for 
equations with random coefficients do not apply, and we therefore use moment techniques 
to study the convergence. 


1 Introduction 

Problems related to wave propagation in random media are encountered in many applications 
that range from imaging the earth’s crust in geophysics m, to communication in underwater 
acoustics [35] or laser beam propagation in the atmosphere mm- The random medium 
often models a complex medium for which only partial information is known. Typically, the 
large-scale variations of the medium (i.e. the background) are known, while the small-scale 
fluctuations (i.e. the heterogeneities) might be too difficult to estimate and are considered 
as random. 

In these applications, waves are generally in a high frequency regime, with frequencies 
sufficiently high so that the interaction of the wave with the fine structures of the medium 
cannot be ignored. From both the theoretical and numerical perspectives, describing the 
cumulative effects of this interaction is a very challenging task. There is therefore a need 
for an approximate, but still accurate, description of the wave propagation. The common 
strategy to attack the problem is based on the high frequency assumption and on asymptotic 
theories of random ODEs or PDEs. There is now a vast literature on this matter, and we 
refer to HU and the references therein for more details. 

In this work, we are interested in two particular asymptotic limits, that we intend to 
perform at once. The first one is the paraxial (parabolic) approximation, which is valid 
when the wave has a privileged direction of propagation and is sufficiently collimated. In 
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the frequency picture, the d —dimensional Helmholtz equation is reduced to the (d — 1)- 
dimensional Schrodinger equation where the time variable plays the role of the variable 
along the axis of propagation. There is a significant gain since a boundary value problem 
is replaced by an evolution problem with lower spatial dimensions. In homogenous media, 
the derivation of the paraxial wave equation is relatively straightforward, and is based on 
asymptotics of the principal symbol of the operator describing the propagation (here that 
of the scalar wave equation). The situation is much more complex when the medium is 
heterogeneous since the interaction with the medium generates some backscattering. One 
has then to resort to some particular features of the medium, for instance small amplitude 
of the fluctuations (5], or oscillatory behavior in nzi, to justify the approximation. 

The second type of limit is of probabilistic nature, and depends on the correlation struc¬ 
ture of the fluctuations. After the high frequency wave has propagated over sufficiently large 
distances in the random medium, it is natural to expect some sort of universal statistical 
behavior to describe the multiple scattering on the wavefield. We are naturally thinking 
here of applications of the (non-)central limit theorem. There is also a vast literature on 
this subject, see for instance mmm- In our context of the paraxial approximation, the 
random medium fluctuations are then asymptotically statistically equivalent to a white noise 
in the main direction of propagation (say z). This holds when the medium has sufficiently 
fast decaying correlations. The limiting model, known as the Ito-Schrodinger equation, is 
studied mathematically in mi- 

When the starting point is the wave equation, or equivalently the Helmholtz equation, 
there are, to the best of our knowledge, only two references on the coupled paraxial-white 
noise limit: in [1], the authors consider the random Helmholtz equation in layered media and 
derive the Ito-Schrodinger equation. Layered media are a nice setting since the dynamics is 
essentially one-dimensional and the transverse variables play little role. In this latter work, 
fluctuations of the medium in the transverse direction are too slow to have a significant 
effect, and the resulting white noise only depends on z. The cumulative effect of the random 
fluctuations on the wave is then a random phase shift driven by a Brownian motion. In EU, 
the medium is much more general, and sufficiently complex to lead to a white noise in z with 
transverse dependence. The cumulative effect is then more complicated and not just a phase 
shift. 

These two references assume that the medium has short-range correlations. It is not 
always the case in practice, as is pointed out in P21123 EH for geophysical problems, wave 
propagation in turbulent atmosphere, or medical imaging. This has then stimulated recent 
mathematical works on wave propagation in random media with long-range dependence [21 
UM 123 ED Eg H3 M- It is shown there that the wave dynamics in such media can be in 
great contrast with that of waves in media with rapidly decaying correlations. For instance, 
anomalous diffusion phenomena were exhibited in nsnsaEu. 

The goal of this paper is to derive rigorously and simultaneously the paraxial and the 
white-noise approximations in the context of random media with slowly decaying correla¬ 
tions in the z direction. Heuristically, the limiting classical white noise is replaced by a 
fractional white noise, leading to the fractional Ito-Schrodinger equation. From the mathe¬ 
matical viewpoint, this is a significantly more difficult problem than the ones addressed in 
in ei. Indeed, in the long-range case, the martingale techniques of mm and standard 
diffusion-approximation theorems for ODEs with random coefficients do not apply. There 
is essentially no general theory in this long-range setting, and we are thus restricted to the 
use of moments techniques which are fairly involved analytically. Note as well that the ex¬ 
istence theory for the fractional Ito-Schrodinger equation is not direct, which leads to some 
additional difficulties in the asymptotic theory. 

Let us be more specific now and introduce the scalar wave equation in the physical space 
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Figure 1: Illustration of the wave propagation model. 


M 3 (the setting could be extended to d > 2, since the techniques used in the paper are 
dimension independent), 


A P- 


1 


c 2 (z,x) * 
equipped with initial conditions 


d?P = V • F 


(t, z, x) G (0, +oo) x M x 


(1) 


P(t = 0, z, x) = dtP{t = 0, z, x) = 0 V(z, x) G M x M 2 . 


Above, the ^-direction will play the role of the main propagation axis, A = d 2 + A x is the 
Laplacian, and A x the Laplacian with respect to the transverse variable x. Here, the forcing 
term F (t, z, x) has the form 


F(t,z,x) := /o(“r - j — )s{z - L s )e z , 

VA 0 tq' 

where 5 is the Dirac measure, e z is the unit vector pointing in the ^-direction, and then 
models a source located in the plane z = L$ < 0, emitting a wave in the z-direction with 
profile /o (see Figure [l]), central wavelength Ao, and transverse width tq. The divergence form 
of the source term is standard in linear acoustics where P would represent the pressure wave, 
see m for instance. Other types of sources could be considered with minor modifications. 
In ([!]), the velocity field is assumed to be given by 

1 ( [l + <rV(f-, f-)] if z6[ 0,L 2 ] 

_ ■= J c o L 1 (c ' c anc l x (z K 2 

c 2 (z,x) | 4 if z € (—oo,0) U (L 2 ,+oo) 

where Co is the background velocity (constant for simplicity), and the random field V(z,x), 
with a stationary covariance, models fluctuations around co in the slab (0 ,L Z ) x M 2 . The 
parameters a and l c represent the amplitude and the correlation length of the fluctuations. 
The main assumption on V is that it satisfies the long-range property in the z-direction, 
which is translated mathematically into a bounded non-integrable autocorrelation function 
which decreases at infinity only as 

E [V(z + s,x)V(s,y)} ~ with fje(0,l), 

z—7-+oo z™ 
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and then 


r+oo 


\E[V(z + s, x)V(s, y)\\dz = +cx>. 


( 2 ) 


We introduce now the scalings, which are similar to these of m where fluctuations with 
rapidly decaying correlations are considered. We assume first that the correlation length l c 
is small compare to the overall distance of propagation in the random medium L z , and we 
denote their ratio by 

l c 

e:=y-Cl. 


Second, we assume that the transverse width ro of the source and the correlation length l c 
are of the same order, 

ro ~ lc- 


This assumption allows for a full interaction of the wave with the transverse fluctuations of 
the medium, leading to a non-trivial transverse behavior. Third, we assume that the central 
wavelength Ao is small compare to L z by taking 


This corresponds to a high frequency regime. With these choices, the Rayleigh length of the 
beam is of order of the propagation distance L z . The Rayleigh length is defined as the distance 
from the beam waist to the place where its cross-section is doubled by diffraction. Hence, 
the beam is still collimated at the exit of the random slab, which is a crucial assumption for 
the validity of the paraxial approximation. In homogeneous media, the Rayleigh length is of 
order Tq/Aq. Therefore, we have for our problem 


Ao 

L z 




Ao 


L, 


'o 

Ll 


This is a parabolic scaling, where the wave oscillations in the z direction are much faster 
than in the transverse direction, which then leads to the paraxial wave equation. From now 
on, we consider the propagation distance L z as our reference scale of order 1, and rescale 
parameters as 

L z = L, Ao = £ 2 lc = £, and ro = e. 

Finally, we consider 


a = e s with s = 2 — $)/2, #£(0,1), 


where S) is related to the decay of the correlation function of V in the variable z as defined 
before. This specific choice of s leads to a nontrivial asymptotic regime in the limit e goes 
to 0. As a result, the wave equation 0 becomes 

AP-i(l + e s v(l)l ( 0 , L) (zj)%P = fo(* -)s'(z - Ls). (3) 

The main result of the paper is the asymptotic description of the pulse front exiting from 
the random section at z = L and around the expected arrival time, which is defined by 

P[{t,x) = P(\ 0 t+ — — —,L z ,r 0 x) = P(e 2 t+ — — —,L,ex). (4) 

Here, the solution is rescaled around the arrival time, and at the transverse scale of the source 
profile. We will show in this work that the process P£ converges in law in C°((0, +oo), L 2 (M 2 )) 
to a process 

pl(t,x ) := J e~ lult ^ UJ (L,x)doj, 
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where satisfies the following fractional Ito-Schrodinger equation 

d^ w (z, x ) = + ik^^z, x)dW H (z, x) = 0, (5) 

where k u := uj/cq is the wavenumber, 

tf w (0,x) := ^e~ iA ^/^f 0 (u;,x), 

with e ~ zL s^x/( 2 k ul ) semigroup of the free Schrodinger equation, and with the convention 

/M = J f{t)e lut dt and f(t) = J f{u)e~ lujt duj. (6) 

Above, Wfj is a fractional field in z with Hurst index 

H := 1 — fj/2 € (1/2,1) 

and will be defined properly further, along with the nature of the stochastic integral. This 
latter integral is of pathwise type, and can be seen as a fractional equivalent to the Ito- 
Stratonovich integral for standard Brownian motions. The function describes the pulse 
deformation, in the paraxial approximation, due to the interaction of the wave with the ran¬ 
dom medium in the section (0, L). The initial condition \Ha;(0) is simply the free propagation 
of the source from z = Lg to z = 0 in the paraxial approximation. In ([5]), backscattering is 
neglected, leading to an initial value problem. As was already observed in different contexts 
in El US for instance, the long-range nature thus leads to a different statistical description 
of the wave than in the classical mixing case of m- In the latter reference, waves are in the 
regime of the central limit theorem, and the resulting Schrodinger equation is driven by a 
standard Brownian field. Here, we are in a different regime where e~^/ 2 Jq V(u/e, x)du con¬ 
verges in law to a fractional field in z with Hurst index H = 1 — fj/2 E (1/2,1). A important 
difficulty in this work is then to justify that a similar type of limit holds for solutions to 
([3]). There are in addition two other main technical points: showing that the coupling with 
the evanescent modes is negligible; these modes exist because of the non-trivial transverse 
frequency content of the random medium; and showing that backscattering can be ignored. 

The paper is organized as follows. In Section |2j we introduce the various assumptions, 
define the stochastic integral and the notion of solution for ([ 5 ]) , and state our main results. In 
Section [3j we give an outline of the proof. The proof is then broken down into the subsequent 
sections. Section [4] concerns the derivation of some important estimates. Section [5] is devoted 
to central technical results about expectation and limits of iterated integrals. Section [6] 
addresses the evanescent modes and Section [7] the backscattering. Section [8] is devoted to 
the convergence to the fractional Ito-Schrodinger equation. Section [9] finalizes the proof of 
the main theorems, and Section |T0| addresses an estimate introduced further. 

Acknowledgment. O. Pinaud acknowledges support from NSF CAREER grant DMS- 
1452349. 

2 Preliminaries and main results 

Throughout this work, we will use the following conventions for the Fourier transform: / 
denotes the Fourier transform w.r.t. the variable t as in ([(jj), and / that w.r.t. t and x, 

f(uj, k) = ( 2 ^\ 3 J f(t,x)e l ( u>t+K ' x ^dtdx with f(t,x) = J f{u,n)e~ l ^ t+K ' x ' > dojdK. 
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2.1 Assumptions 

The source term. We suppose that fo(u>, k) is a bounded function with compact support 
in both variables, and even in the variable c o. We assume moreover that it is supported away 
from zero w.r.t. oj, that is there exists io c > 0 such that 

(-w c , w c ) n supp u f 0 (w, k) = 0, Vk. g R 2 . (7) 

The latter assumption essentially means that the source is shorthand. Larger bandwidths 
could be included by direct modifications of the proofs. 


The random field. We construct the random field on a probability space (H,T, P) and 
in the Fourier space as follows: the field V is the Fourier transform of a random measure 
V(z, dq), i.e. 


V(z,x) = 


0 -iq-x 


V( z , dq). 


We define V(z, dq) sufficiently explicitly in order to be able to carry on the calculations. Let 
then S C R 2 be a bounded domain, symmetric around the origin (S = —S), included in a 
ball B(0,rs)- The domain S will be the support of V since we want the largest transverse 
frequency to be at most of order e _1 (after rescaling). Let also 3% be a real-valued mean-zero 
Gaussian random field on [0, +oo) x S, continuous and stationary with respect to the variable 
z, and such that B^(z, q) = 3^(z, —q). Its covariance function is given by 


K[Bsj(z + zo,qi)Bsj(zo,q 2 )] ■= r^{z)R{q 1 ,q 2 ), 


where R is assumed to be a continuous positive symmetric and bounded function such that 


0 < R{qi,q 2 ) < R(q, q) = 1, V(g, qi, q 2 ) G S X S x S. 


Besides, r ^ is a continuous even function bounded by r^(0) = 1 and 

rsj(z) ~ -j- with Sj G (0,1). (8) 

z—>-h oo Z 3 * 

Hence, is not integrable at the infinity. Let then © be a smooth odd function satisfying 
for all Z £ N, 

sup |0^( , u)| < Cq, (9) 

uSK 

where Qd) stands for the Z-th derivative of 0, and consider @(Bsj(z, q)), which is not a 
Gaussian variable. Introducing a random measure m(dq), supported on S, independent of 
the random field 3%, and whose properties are defined below, we write V as V(z,dq) := 
m(dq)@(Bsj(z,q)), so that 


V(z,x) = [ m(dq)e iq ' x @(B^(z, q)). (10) 

Js 

We suppose that m*(dq) = m(—dq), with bounded associated total variation measure \m\, 
that is, almost surely, 

\m\(S)<C m , (11) 

for some deterministic constant C m > 0. This yields in particular that V is real and bounded, 
and therefore that the velocity field cannot take negative values for e sufficiently small. We 
suppose moreover that m is stationary, 


E[m(^i)m(</? 2 )] = / m{dq)(p 1 {q)^ 2 {q), 
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Figure 2: Examples of realizations of V. Here, V is obtained via (10) with 0(x) = sin(lOx) 
and S = [—15,15]. The measure m is as in (12), where the Uj are uniform in [—1,1], the qj 
are chosen with a discrete uniform distribution among the points of a uniform discretization 
of S, and aj = 1 for j < 100 and zero otherwise. The field B$ is obtained via a similar formula 
as (16) where the e n are cosines, the j3 n behave like n -2 and Wh,u is replaced by a Gaussian 
process with autocorrelation r^{z) defined as the Fourier transform of l(-i5 J i5)(fc)/l^| 1- ‘ 
From left to right, Sj = 0.9, 0.5, 0.1. Observe the arising of long-range correlations in the 2 
direction as Sj decreases, that is as r$ decreases slower at the infinity. 


where m is a positive measure on S with finite mass and the (p t are smooth functions. This 
yields that V has a stationary covariance in both z and x, which is a common assumption in 
applications. An example of such a measure is the following: 

m(dq) = aj(UjSqj + UjS- qj ), ( 12 ) 

j> o 

where (aj)j>o E Z 1 (N,M) is deterministic, (Uj)j> o E C N and (qj)j> o E are independent iid 
sequences of random variables with appropriate distributions, and the Uj have a zero mean. 

Note that E\V{z, x)] = 0 by symmetry, and in the same spirit as [ 251 Lemma 1], we show 
in Proposition |5.1| of Section [5] that V itself satisfies the long-range property 

E[V(z + zo,x)V(z Q ,y)] ^ ~ R 0 (x - y ) with := — ( / Q(u)ue~ u/2 duj , 

Z / i CXD £ ZjIT ' J —oo ' 

and 

Rq(x) := / m(dq)R(q, q)e~ iq ' x = f m (dq)e~ iq ' x . (13) 

Js Js 

This implies that ([2]) is satisfied. Examples of realizations of V are given in Figure [2j 

The limiting field Wh of ([5]) is heuristically obtained as follows: the scalings in V and 
the long-range behavior act in a such a way that only the linear part in 0 is not negligible, 
and such that Jq q)du is well approximated (in distribution) by a fractional 

Brownian field in z. Hence, Wh is a random field with covariance operator given by 

E[W H (z 1 ,x)W H (z 2 ,y)] = ^ ^ ( z\ H + z\ H - \ Zl - z 2 \ 2H )R 0 (x - y ) (14) 

for all (zi,z 2 ) E [0,+oo) x [0,+oo) and (x,y) E M 2 x M 2 . The construction of Wh and the 
definition of the stochastic integral are given in the next section. 


2.2 Stochastic integral and fractional Ito-Schrodinger equation 


The stochastic integral with respect to a fractional Brownian motion obtained here in the 
limiting process is of pathwise type, and is defined according to the work of Zahle |36| . We 
start this section with the construction of the fractional field Wh with covariance operator 
(14), which is used thereafter to define the stochastic integral. Finally, we give the definition 
of a solution of ([5]) before stating the main results of the paper. 
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Fractional field. A one dimensional standard fractional Brownian motion with Hurst index 
H, on a probability space (fi, T, P) is a centered Gaussian process bn with covariance 

E [6# (it) &#(?;)] = -(u 2 ^ + v 2H — \u — v\ 2H ), V(«, v) £ [0, +oo) x [0, +oo). 

Moreover, such a process admits the following spectral representation 

n Aru _ i 

b H (u) = C]j 2 J . r ^ H _ 1/2 w(dr), (15) 

with Ch = HT(2H)sm(TrH)/ir, and where w{dr) is a complex Gaussian random measure 
such that w*(dr ) = w(—dr) and 


E[rc(dr)u)*(ds)] = 5(r — s)drds. 


The construction of the fractional field with covariance operator (14) is done in the Fourier 
domain. Let (H / A,ra)n>i be a sequence of independent standard fractional Brownian motions 
on the probability space (fl, T, P) with Hurst index H. Using the fact that 


Q : L 2 (S) 



L 2 (Sl 

fs dpR(p, q)tp(q) 


is a positive self-adjoint trace class operator [Hj Corollary 4.4], their exist a sequence (e n ) n >i 
of orthonormal eigenvectors and a sequence (/3 n ) n > l of positive eigenvalues for Q. Therefore, 
the Gaussian random field 


B H (z,q) := ^2 VJ3n e n(q)W H , n ( z ) 

n> 1 


(16) 


defines a infinite-dimensional standard fractional Brownian motion on L 2 (S), and then 


Wh(z,x):=cth [ m(dq)e iq ' x Bii^z^q) with 

Js 


(j 2 H = 


C f 


H(2H — 1) 


(17) 


defines a random mixture of fractional Brownian fields with Hurst index H and covariance 


operator given by (14). It is interesting to note that our asymptotic noise model is not 


Gaussian. In fact, for medium perturbations with mixing properties the asymptotic noise 
is always Gaussian m na E|. However, for medium perturbations with slowly decaying 
correlations, it is not necessarily the case. In a one-dinrensional wave propagation setting, 
it has been observed in |26] that the asymptotic noise model is not necessarily Gaussian 
if the initial medium fluctuations have non-Gaussian statistics. In our context, the non- 
Gaussianity of Wu is due to the random measure m, which has been introduced for the 
covariance stationarity of V. Nevertheless, the form of Wh allows the use of the Gaussian 
properties, and then the use of |36j to define the stochastic integral in (|5]). 


Stochastic integral. We follow here the approach of 
space 

IU Q (0, L, B) := j?/' G C°([0, L\, B) such that HV’lla,# < +ooj, 


. Let us consider the Banach 


| a,B ■= sup 
ze[o,L] 


* B + 


z) - (/>(u)\\ 


B 


— 7/l Q +l 


{z - u) 


du 


with 










and where B is a given Banach space. Moreover, for (3 E (0,1), let us denote by M([0, L\, B ) 
the Banach space of f3 -Holder functions on [0, L\ with values in B , equipped with 


\\<t>\\p,C,B ■= sup \\4>(z)\\ b + sup 

ze[0,L] 0 <v<u<L 


ll</M - 4>[v)\\b 

(u — v)P 


One can remark that for a E (0,1/2) 


W\a,B < L^Ulh.^B SO that C 1 ~ a ([0,L],B)cW a (0,L,B). 


(18) 


Now, for a E (0,1) and z E (0, L), let us introduce, for a real-valued function /, the so-called 
Weyl’s derivative given by 


:= 

D a L-f(z ) := 


1 


r(l-a) 

(-1) Q 


m 


+ a 


f{z) - f(u ) 


— 7/1°+1 


du 


r(l - a) L (L-z) 


f(z) 


+ a 


(■ z-u ) 

rL f(z) - f{u) 

z (u - z) a+1 


du 


whenever these quantities are well-defined, and where T(u) = r u ~ 1 e~ r dr is the gamma 
function. Following [3S], the generalized Stieljes integral of a function / E C^([0,L],M) with 
respect to g E C^([0, L\, M), with v + n > 1, v > a, and fj, > 1 — a is defined by 

[ L fdg := (-1)“ [ L DfrfMD^gL-Mdu, (19) 

Jo Jo 

where gL-(u) '■= d(u) — g{L~). The definition does not depend on a, and we have 

Jo ^ d9 = Jo 

Moreover, according to |28], this integral can be extended to more general classes of functions 
thanks to the relation 


where 


and 


A a(g) ■= 


fdg 


I/Ml 


< Wf\\a,iK(g), 

I/M - /Ml 


+ 


dv ) du, 


I o (u - u) a+1 


sup \D t _ a g t -(u)\. 


r(i Qi ) r ( oi ) o <u<z<l 

Consequently, this integral is well-defined as soon as / E W“(0, L,M) and A a (g) < +oo. 

As a result, for a random function F E W a (0, L, L 2 (M. 2 )), the stochastic integral with 
respect to the fractional field Wh, 

[ F(u,x)dWH{u,x), 


is defined by (19) almost everywhere in x and P-almost surely. In fact, we have for a E 

< C\\F\\w a (0,L,L 2 (R 2 )) SU P Aa/IFffM), 


(1 — FJ, 1/2) and for all z E [0, L\ 
F(u)dW H {u 


L 2 (R 2 ) 


with 


E 


sup A a (W H {x)) < / E[\m\(dq)]E[A Q (B H (q))} < C m supE[A a (B H (q))} < oo, 
xeiR 2 J Js ges 


as will be proved later in Lemma 8.1 
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Fractional Ito-Schrodinger equation. The notion of solution for the Ito-Schrodinger 
equation ([ 5 ]) is made precise in the following definition. First, let us introduce some additional 
notations. Let k G N*, and let us denote by H k (M?) the fc-th Sobolev space on M 2 . Consider 
moreover Wjf(0,L) := W a (0, L, H k (M. 2 )), equipped with the norm || • ||j an d the 
complete metric space 

W“(0,L):= fl W£(0,L), 

k£ N* 

equipped with 

1 


dn 


VV>,V0 :=E^( 1A 


I a,H k (R 2 ) 


k> 1 


Definition 2.1. Let H G (1/2,1) , a G (1 — H, 1/2), and Wh be the fractional field defined 
by ©■ We say that G Wff(Q, L) is a pathwise solution of © if, with probability one, 
for all (z,x) G [0, L\ x M 2 , we have 

i f z f z 

^ u (z,x) = ^ L0 ( 0,x) + —— / A x ^ 0J (u,x)du +ik u ^ UJ {u,x)dW H {u,x). 

2 k w Jo Jo 

In other words, a solution of ([5]) is a pointwise solution of this equation for almost all 
realizations of the randomness. We will see later that a solution to © has automatically 
Holder regularity 


^GC(0,L):= fl C 

fceN : 


H-e 


([0 ,L\,H k (R 2 )) with 6 = H + a- 1. 


Here, Cffi y (0, L) is a complete metric space equipped with 


dH-0,C,oo(4>, VO EZ fik A HV 1 “ V , ll/f-e,C,H fc (R 2 )) 


k> 1 


so that C^f~ e {f),L) C W^(0,L) according to (18). The solutions we define here are classical 
solutions in the standard terminology. It is not completely trivial to construct less regular 
solutions to ©, which is required for non-linear problems, see [213]. Indeed, the standard 
technique is to use the mild formulation, and then treat a term of the form 


S(z — u)^t u (u, x)dWn{u, x), 


where S is the Schrodinger semigroup. As explained before, some Holder regularity in u 
is needed in order to make sense of the integral. Since the semigroup is not sufficiently 
regularizing, this regularity in u has to be exchanged for some regularity in x on 'L a; , and 
the fixed point procedure cannot be closed. This is not a problem in our linear setting where 
we can iterate the stochastic integrals and suppose that the initial condition is C°° in x. A 
different strategy has to be adopted in the non-linear case |29|. 

Note that the stochastic integral here is the fractional equivalent to the Ito-Stratonovich 
integral for standard Brownian motions, and as such satisfies the classical integration by 
parts formula. This then formally yields the conservation relation, for all z G [0,L], 


1 , 


I^u;(^)||l 2 (M 2 ) — II ^a)(0) I|l 2 (R 2 ) — xII/o(^)||l 2 (R 2 )- 
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2.3 Main results 


We will actually not work directly with the process Pf given by Q, but rather with an 
approximate process p e L defined by 


p e L (t,x) := p(e 2 t + -—— ,L,ex), 
v c 0 > 


( 20 ) 


where p solves the wave equation © with artificial absorption a £ , 
1 


A P~ ~a ( 1 +£SV {~' t) 1 ^ ,L){z))d?p + ia £ p = /o(t 2 , ^) s '( z ~ L s ) 


and vanishing initial conditions. Thanks to the estimate below, proved in Section 10 

C T aJ 2 

sup ||P(t, £■) - p(t, z, £■) || L 2 (R 2 ) < - VT > 0, 

(t,z)e(0,T)xR e 


( 21 ) 


it is equivalent, from the viewpoint of convergence in law (see [BJ Theorem 3.1 pp. 27]), to 
consider p e L instead of Pf by choosing a £ = o(e 2 ). The main theorem will be hence stated in 
terms of p e L . The introduction of p is an important point since the absorption term provides 
us with straightforward estimates in L 2 ((0, L) xl 2 ), that would require much more work with 
the process P. These estimates are not uniform, but sufficiently tamed, and are exploited 
throughout the paper. 

We will mostly work in the frequency domain, and in order to take Fourier transforms in 
time, we extend p to negative times by setting p(—t, z, x) = p(t, z, x), for all t > 0. 

The main result of this paper is the following theorem, that states the convergence of the 


pulse (20). 


Theorem 2.1 (Convergence result). The family (pf,) £ g(o,i); defined by (20), converges in 
law in the space C°((— oo, +oo), L 2 (M 2 )) D L 2 ((—oo,+oo) X M 2 ) to a limit given by 


p° L (t,x) = j e tujt ^ u (L,x)du, 

where 'Fa, is the unique pathwise solution to the fractional Ito-Schrodinger equation ([5]) . 

The second theorem below is a by-product of the proof of the main theorem, and provides 
us with some interesting properties of the solutions of the fractional Ito-Schrodinger equation: 
existence and uniqueness, conservation of the energy, approximation by a smooth process 
which can be expanded in terms of scattering events, and approximate formulas for moments 
of any order. 


Theorem 2.2 (Properties of the fractional Ito-Schrodinger equation). We have the three 
following statements: 

1. The fractional Ito-Schrodinger equation © admits a unique pathwise solution for 
all a € (1 — H, 1/2), which satisfies 

II^(^)||l2(r2) = -||/o(^)||l 2 (r2) Vz €E [0,L]. (22) 

Moreover, € C^fi e (f),L) for all 6 6 (0, H — 1/2). 
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2. For all 6 E (0, H — 1/2), the process can be approximated by 

^u(z)= lim V£(z) with v£(z) = /o(w, ■) + Y, 

^+°° 

where the limit holds in C^f~ e ( 0, L) m probability. Here, we have in the Fourier domain 
*Z' A (z, k) := (ik u ) n f m((iq (n) ) [ 


n 

m=l 


fA e ir m u 


I—A \r n 


, H _ 1/2 w(dr m ,q m ), 


where u (n) := (ui,.. .,u n ), q (n) := (qi, ...,q n ), m(dq (n) ) = m(dq 1 ) . ..m(dq n ), Q (n) : = 
(Qo, • • ■, Q n ), with Q m := k - qi - q m , and 

1 n 

G n ( U W,QW) := — £ (|Qm-l| 2 - \Qm\ 2 )u m - 

ZiHji 


171= 1 


Moreover, 

A„(z) := {(iti, • • • ,u n ) e [0,z] n , s.i. 0 < Uj < Uj-i Mj G {2,..., n}}, 
and (w(dr, q)) q£ s is a family of complex Gaussian random measure defined by 


w 


(dr, q) := 


lp n CsjT(2H- l)sin(Trtf) 


n> 1 


7r 


e n (q)w n (dr), 


(23) 


where (w n (dr )) n >i *s i/ie family of independent complex Gaussian random measure in 
the spectral representation (15) of the family {W n ,H)n> l introduced in (fT7|) . 


5. We Ziaue 


E 


Mi 


M 2 


nw*) n 


j 1=1 


J2=l 


= lim E 

A —>+oo 


Mi 


i 2( R 2(Mi+M 2 )) 

M 2 


n < Jt a) n 


,11' / J.J. ^2,j 2 ' 7 ' ' / I,2( R 2(M 1 +M 2 )J ’ 

3 1 = 1 J2=l 


/or a/I (Mi, M 2 ) 6 N 2 , frequencies (w i)3 -)(i ! j)e{i,2}x{i,...,M i } and p 6 L 2 (R 2 ( Ml+M ^). 


The conservation relation (22) is a consequence of a negligible backscattering and shows 


that the energy of the pulse is conserved at the end of the random section, that is 

IIPlII-L 2 ((— oo,+oo)xIR 2 ) 7) ll/oII l 2 ((— oo,+oo)xR 2 ) ■ 

Note also that the convergence in the second point holds in W^(0,L) with a = 1 — H + 6 
according to (18), and that becomes smoother in z as H increases, which is expected 


since the regularity of the fractional brownian motion improves with H. Moreover, we will 
see further that the process ’ll/) is the solution to a fractional Ito-Schrodinger equation with 
a regularized fractional white noise, and as such enjoys some regularity properties (w.r.t. z) 
that are convenient in justifying formal computations, in particular the calculation of the 
moments as in item (3) above. Moments are important for instance in imaging applications, 
where they help quantify the stability of reconstructions with respect to changes in the 
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random medium, see e.g. mmmm- The series expansion in item (2) is the classical Born 
approximation, see e.g. [2H Section 17.2], 

Note finally that \k a; satisfies various formulations of ([5]), for instance 


$u(z, «) = ^u>(0,«) - 


l\K\ 

~2ku 


yiu(u,K)du + ikuGH / m(dq ) / $ u (u, k - q)dB H (u, q) 


in the Fourier domain, or the mild formulation 


$ w (z,k) =e~ i ^ 2z ^ 2k ^4f LJ (0,K)+ik UJ a H f m(dq) [* e -i|«l a (*-«)/( 2 fct-)^ w ( u , K - q )dB H (u, q)du, 

Js Jo 

where the relation between Wh and Bh is given by ©■ 

The rest of the paper is dedicated to the proofs of the theorems. We will focus on Theorem 


2.1 since its proof contains that of Theorem12.2 1 We begin with an outline in the next section. 


3 Outline of the proof 


The proof starts with recasting the wave equation in the Fourier domain in all variables, 
except z. Defining for this the Fourier transform in time 

rM = / g(t)e iuJt/£2 dt with g(t) = j F^e-^du, 

that accounts for the high frequencies generated by the source, we obtain from © the 
Helmholtz equation 


+ A x p £ u (z,x) + ^(l+£ s V^, j)l ( 0 , L / e )(z)')pZ,(z,x) + ia e pfj(z,x) 
= fo(u, -)s\z - L s ). 


(24) 


We will construct solutions to (24) in section [4j and show that satisfies the required 
regularity to justify all the calculations. Taking the Fourier transform w.r.t to x of the 
wavefield rescaled around the propagation axis pl,(z,£x), that is 

pH z ,k) := 77^2 [ p £ UJ {z,£x)e lKX dx, 


(2tt) s 


we find 


+ K ) + ia efui z > K ) 


£ 2 M 2 ' 


+ £ s A kl f m(dq)v(-,q]p e w (z,K-q) = fo(uj,K)S'(z - L s ), 

J s & 


(25) 


where k u = ui/cq is the wavenumber, and V = @(H^) according to (10). Following the 
standard terminology in absence of absorption (a £ = 0 ), when the wavevector k satisfies 
|k| < \k u \/e we will refer to the corresponding mode as a propagating mode. These modes 
can propagate over large distances. When the wavevector satisfies |/c| > \k u \/£, we refer to an 
evanescent mode. The proof is then based on a decomposition of pfj into right- and left-going 
propagating modes, see Figure [ij with amplitudes af, and 5^, respectively. These amplitudes 
are defined as the solutions to 


F u (z,k) = -p^=(o5(2,«)e <fc “ Ae '"W*/ ea +S5(«,K)e- ifc “ Ae '"W z / e2 ) 

\J ^e,u{ k ) 


(26) 


dzPlj( z ,K) = 


tkuy/KM K yk„\ e , u ( K )z/e* _ , (27) 
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Figure 3: Illustration of the right- and left going mode amplitudes. 


for |/c| < \k u \/e, and where 

:= £ 2 \k,\ 2 /kl + ia WtS , with 


OL^S 

auJ ’ £ = ~T2~- 

' h t,) 


(28) 


We consider here the principal square root for complex numbers, namely the square root 


with positive imaginary part. Using (|25|) we find the following coupled-mode equations 

d z 


Oy(z,K) 

K{z,k) 


+ £' 


{|«— q\<\kuj\/£}r\S 

i/krj 


m (dq)H^(z, k, q) 


dy{z,K- q) 
b e u {z,K- q) 


s -2 


2 \JK,gj(k) Al«-?l>M/e}nS 


m(dq)V(z/e,q ) 


p—ikujXz^i^z/e 2 
_ p ik w X S}U (K)z/e 2 


Pw(z,K- q), 
(29) 


for |ac| < \k u \/e, with 

H„{z,K,q) = e s ~ 2 


ik u 


--V(z/e,q) 


2-y/ ^£,u(q) 

e ik u (X e ,ut(q)—X Sl u(K))z /£ 2 e —iku{X et u(q)+X e ,u(n))z/E 2 
_ p iku(Xe,u(q)+Xe,u(n))z/E 2 _ p -ik UJ (X S}LLj (q)-Xe,ui(ii))z /£ 2 


(30) 


The system (29) is equipped with the boundary conditions 


a^(0,K) = 




e - ik “ x *’“ Ls / £2 f 0 (uj,K) and SP(L,/c) = 0, 


(31) 


where the first condition represents the (known) amplitude of the wave coming from the left- 
homogeneous half-space and entering the slab (0 ,L), and the second condition implements 
the fact that no wave is entering the right-hand side of the slab. These conditions will be 
investigated in more details in Section [4] Despite its formulation, the system (29) is not an 
initial value problem (IVP), but rather a boundary value problem. The limiting problem will 
nevertheless be shown to be an IVP. 
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Note that in the homogeneous case where V = 0, the justification of the paraxial approx¬ 
imation is straightforward. Indeed, in that case afj(z, k) is constant for z > L$ and therefore 
equal to a^(0, k) defined in (31), so that the right-going wave a is given by 

z- L s 


p e Jt,x ) =p(e 2 t-1 - -,z,£x] 

V Cn / 


Co 

—iut —ii tx e iku, (.\ e 


'W-l ){Z - L s)/ £2 f 0 {uj, K ) d LU d K, 


and therefore, pointwise in (t, x ) 


lim pF z (t,x) 
£—^0 


,-iwt B -inx-i$-(z-Ls) 


fo(u], K)diodhi 


= - e 


—iu)t 


-{z L S )A x j^, .)] (x)du, 


where e * 2A z/( 2fe A j g ^he semigroup of the free Schrodinger equation. 

When V is random, the core of the proof is an asymptotic analysis of the amplitudes a £ u 
and b^j solutions to (29). There are several steps, most of which involve computing moments 
of the form 


E 


[n 

3 =1 


(32) 


The random variables Q(B^(u/e, q )) are not Gaussian, but we will see in Section [5] that they 
asymptotically behave as Gaussian variables. This is based on the ideas of (34] : suppose that 
n is even; the even function 0 is then projected on the basis of the Hermite polynomials, 


which gives rise to a series of terms in (32). The leading one corresponds to the product of 


the first order terms in the Hermite expansions, and consists of the product of n/2 moments 
of order 2 of B Any other term involves at least n/2 + 1 moments of order 2 which, after 
integration, the use of the scaling Uj/e and of the long-range depen dance property, leads to 
negligible contributions. The fact that a fractional Brownian field is obtained at the limit is 
a direct consequence of the asymptotic behavior of the correlation function (J8l) . 


Owing the technical results on terms of the form (32), we can then proceed to the analysis 
of a £ u and b £ u . The first step is to show that the coupling between propagating and evanescent 
modes can be neglected. This is actually a fairly subtle point. The fact that evanescent 
modes are expected to decrease exponentially as z increases in (0, L) cannot be exploited 
close to the transition propagating/evanescent modes. Indeed, in this case, X £tU is too small 
around the transition, which essentially yields exponentials in (26) of order one. Our method 
then goes as follows: we start from (29) with |/c| < \k e \/e, and therefore only consider 


the propagating modes; the first term in the r.h.s corresponds to propagating modes with 
momentum k — p scattered to propagating modes with momentum k after interaction with 
the random medium; the second term corresponds to evanescent modes with momentum 
k — p scattered to propagating modes with momentum n. We will prove the convergence of 
(d/,, bfj) in the distribution sense, which will limit |k| to some bounded domain independent 
of e, say |ac| < R. In this case, the second term is zero, since only propagating modes 
with large wavenumbers of order \k u \/£ (which is of order e _1 since k^ is bounded from 
below independently of e according to assumption (|7|) are coupled to the evanescent modes. 
This naturally does not mean that evanescent modes have no influence on the propagating 
modes, the coupling appears in the first term of the r.h.s via modes with larger and larger 
wavenumbers that get closer to the transition. 

As an approximation, it is therefore natural to introduce the following system, that only 
describes propagating modes, for all |«| < \k^\/e, 


AUz,k) 

B £ u (z,k) 


a £ u ( 0 ,k) 

k,( 0, K ) 


+ 


du 


'{|K-q|<|fc„|/£}ns 


m(dq)H £ (u,K,q) 


A £ u(u,k- q) 
Bl{u,K-q) 


(33) 
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extended to 0 for |«| > \k u \/e, and where 


H^{z,K,q) = 


IE 


s—2 




V(z/e,k- q) 


(AJ (<j) — AJ w (k))z/ e 2 


(q)+\ r e (K))z/e 2 


Jk u {\ r e ^{q)+\ r e ^{K))z/e 2 e -ik il ;{\ r E ^{q)-\l^{K))z /£ 2 


(34) 


with 


K,M := \A - e 2 l«l 2 /e 


The important fact is that evanescent modes are absent in the system above, a minor point 
is that H J is approximated by H^, where the absorption a e is set to zero and the X EjLJ in the 
denominator in are set to one. We want to show that (A^, Bff) is a good approximation 
of (&£,,&£,). For the sake of clarity, let us assume that H is replaced by in (33) since 


this is not the main issue. Integrating (29) in z and taking the difference with (33), we then 
need to prove that two homogeneous solutions to (33) for |/c| < R, that is solutions with 
the first term in the r.h.s set to zero, are close to each other as e -> 0. The solutions are 
not equal for e fixed, otherwise there would not be any evanescent modes in (29), which 
is obviously wrong. The main difficulty is that the system is not closed: even if we only 
consider modes with |ac| < R in the l.h.s., modes with larger wavenumbers are involved in 
the r.h.s.. We will use the following observation to overcome this issue: two modes with very 
different wavenumbers, say one with \k\ < R and the other with \k'\ S> R, are related to each 
other only after a large number of interactions with the medium; since there is some loss of 
amplitude at each interaction (due to scattering and not the artificial absorption a e ), the 
coupling between the |k| and the |ac 7 | modes is expected to be small. This idea can be seen as 
a form of asymptotic closure of the non-closed system, and is implemented in the following 
proposition, proved in Section |6j 

Proposition 3.1 (Coupling with evanescent modes is negligible). For all z £ [0,L], for all 
p > 0, and for all test function 0 £ Cq°(M 2 ) x Co°(M 2 ), we have 


£—^0 


(/ 

<%,(*) 


'K(z) 

l \ 

K{z) 


K{z) 


L 2 (R 2 )xL 2 (R 2 ) 


> M = 0. 


An easy consequence of this result is that for all p > 0, T > 0, and test function cj) with 
£ Co°(M 2 ), we have 


limP( sup \(p e L (t) -p\ tL (t), > fi) = 0, 


(35) 


where p e L is given by Q and 


— iijjt -IK-X ' 


p £ 1L (t,x) := / / dudne lut e 


,-ik^lL-Ls)tE 2 


A £iW (k) 


i^(L,K)e a '“ As '“ (K)L/£2 


Let us remark that the left-going mode amplitude is not involved in p\ L because of the 
boundary condition 6^,(L) = 0 (no wave is coming from the right homogeneous space). Since 
pf—PiL converges in probability to 0, it is then enough to investigate the limit in law of 
Pi l prove Theorem 2.1 (see P, Theorem 3.1 pp. 27]). Even if this latter convergence only 


holds in a weak sense, energy estimates given in Section [4] will allow us to obtain the strong 
convergence in L 2 (M 2 ). 

The second step of the proof is to study p\ L and therefore the couple The 


system (33) is closed, but the backscattered mode amplitude 6^(0) is unknown and cannot 
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be considered as an initial condition. This issue is corrected by introducing the propagator 
P®, defined as the solution to, for all |k| < |/c w |/e, 


^d z ,^ r ) = 


S(k — r) 0 
0 5(k — r ) 


r z r 

+ / 

du / 

Jo 

HI 


{|K-g|<IMMruS 


m(dq)Hu{u, k, q)Vu{u, k - q,r), 


where the initial condition is the identity operator. The propagator is extended to zero for 
|k| > \kuj\fe. According to symmetries in H^, the propagator can be decomposed as follows 




V% £ (z,K,r) Vu' £ (z,K,r) 
V b f{z,K,r) V% £ (z,K,r) 


where (' P^ £ ,V^ £ ) is the solution of 


V% £ {u,k- q,r) 
V% £ (u,k- q,r) 

(36) 

The term V% ,£ describes scattering to the same direction of propagation, while V b,£ describes 
scattering to the opposite direction. We then find the following relation between the right 
and left going modes in terms of the propagator: 


V^ e {z,K,r) 

V b j £ {z,K,r) 


S(k — r) 
0 


du 


'{\K-q\<\k u \/e}nS 


m {dq)H^(u, k, q) 


A £ u (z,k) 

K{z,k) 


= J K( z ’ K > r ) 


Ow(0 ,r) 
6^(0, r) 


dr, 


Mz £ [0, L\. 


The expression of the wave exiting the random section is thus, after integration against a 
test function £ Cq°(M 2 ), 


{p\,lA)A) -=P £ a,dtA) +Pb,L{t,(f>), 


with 


P^l^A) ■= jj dudne luJ AdAAd j drV a f{L, re, r)<j> e u (r) 

Pb^A) ■= jj dLjdre~ ,wt e tkuLs/£2 bl J (0,r) j dnVt’ £ {L, k, r)Ad K AA)- 
Here, we have defined the following functions, 


Add : = 


\fKJA) 


f 0 (u,r)e- ik ^ x ^-V L s/ £ \ Add ■■= 


Jkui (A e , u (k)—1 )L/e 2 


A £)W (fc) 


Before describing the asymptotic behavior of the propagator, we need to introduce a few 
more notations: 


:= J V A( L ^,dAd r )dr, V*’^l(L,r) := J Vd{L, k, r)Add<f>A)dn 

Add ■= fod,r)e %M2L s/^_ 


(37) 


The next result shows that the backscattering is negligible and provides us with the leading 
term in the propagator. 

Proposition 3.2 (Backscattering is negligible). We have the following two statements: 
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1. For all r] > 0, we have 


So P (/ K^,( l )\\l 3 (sl 3 )^ >v)=0, 

where V*J is extended by 0 for |/c| > Ifc^l/e. 

2. For all r) > 0, we have 

limP( J K%(L) - Xf(L)\\ L 2 {R2) du > v) = 0, 

where Xf is the solution to 


d z X*{z,n) = f 

■F 


{|re—g|<|fc w |/e}n<S 


,i(z, k, q)Xf(z, n - q) with Xf(0) = <j>\ 




for 

/-ryCLj£ 


k| < \k w \/e, and H J is defined by ( |34| ) l 

and Xf are extended by 0 for |k| > \kj\/e. 


(38) 

is the (1,1) entry of Hf). Ftere, 


Proposition |3.2| is proved in Section [7] The proof is based on a series expansion of the 
propagator and on the fact that the coupling between right and left going modes appears via 
an oscillatory integral. Note also that neglecting the backscattering leads to an IVP on Xf. 

The last step is to characterize the limit of Xf. With the same arguments as before, we 
only need to investigate the convergence in law of Xf to prove Theorem 2.1 The proof, given 


in Section [8j is also based on a series expansion of Xf and on the computation of the limiting 
moments of Xf. We will need the following functional spaces: for k £ N*, let us denote by 


ttk := j</> = $ with if G iL fc (M 2 ) j with 


|2 

\H k 


/(1 + \K,\ 2 ) k / 2 \(f(K)\ 2 dn. 


Consider also Wjf(0,L) := W a (0, L,FLk), equipped with the norm || • || Q) fc := || • \\ a ,H k i an d 
the complete metric space 

W*(0,L):= f) Wg(0,L) equippedwith d a ,oo{F, VO ■= ^ i (1 A \\4> ~ V’lkfc) • 
feeN* k>1 

We finally introduce the complete metric space of 'Hfc-valued functions with Holder regularity 

Cg~°(0,L):= H C H ~ e ([0,L\,H k ), 

fee N* 


equipped with 


dH- 0 ,c,oo((t>, VO ^ A ^ ~~ ^\\H-efi,u k ) 


k> l 


We will also use the spaces C^(0,L) and C“(0,L) with immediate definitions. Let us recall 
that according to (18), we have C^~ e (0,L) C Wff(0,L) for 9 = H — a — 1. The convergence 
result is the following: 


Proposition 3.3 (Convergence to the fractional Ito-Schrodinger equation). We have the 
three statements below: 


1. For all M 6 N* and frequencies (wi,...,wm), the family (X^^L),..., X% (L)) con¬ 
verges in law in L 2 (M 2M ) to (X Ul (L),... ,X Um (L)). Ftere, X u is the unique pathwise 
solution in L), for all a G (1 — H, 1/2), of 

X u (z, k) = (k) + ik u OH [ m(dq) [ due~d\ K - q \'-\F 2 ) u /{^) K _ q )dB H (u, q ), 

Js Jo 

(39) 
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where Bn is defined by (16) and an by & Here, the stochastic integral is defined P- 
almost surely pointwise in n and q. Moreover, X u G C^~ 6 (0, L) for all 0 G (0, H — 1/2), 
and 

11^(20 IIl 2 (r 2 ) = 11^(0) ||l 2 (m 2 ) = 2 H/oC^’ ')lli 2 (K 2 )- (40) 

2. We have for all (Mi, M 2 ) G N 2 , all frequencies (wij)(jj)e{i,2}x{i,...,Mi}> an d test function 
ip G L 2 (R 2(Mi+M2 )), 


lim E 
£—^0 


Mi M 2 

TT xl (l) TT XL (l),<p, 

11 ^1,31 v ' 11 a; 2j 2 V n “/j2(|2(M 1 +AJ 2 )) 

jl = l 32=1 y ; 

A/ 1 A /2 

ii=i j2=i 


t 2 (t 2 («l+M 2 )). 


5. Therefore, the process defined in the Fourier domain by 

*Uz,k) :=e -^lV(2 '*->)x u (z,k) 

satisfies all the requirements of Theorem 2.2 with Wh defined by ([ 13 - 

In order to identify the moments of X u with the limits of those of Xf, we will identify 
the moments of each term in a Duhamel expansion of X :jJ . The main technical difficulty 
is to handle the fractional stochastic integral and to justify the calculations, in particular 
exchange of expectation, limit, and integration. We will for this proceed by regularization, 
and start by constructing an approximate solution X A that solves (39) with Bh replaced by 
Btf defined by, for q G S, 


B H(u,q) := C]j 2 ^ V %.e n (q ) 


n> 1 


r A e iru - 1 
_A irlrl -^ -1 / 2 


w n (dr), 


with Ch = HT(2H) siifiirH) /tt, and where (w n (dr)) n > 1 is the family of independent complex 


Gaussian random measure given in the spectral representation (15). Since Bfi is C°° w.r.t. 


to u, the integral in (39) is now simply a Lebesgue integral, and computations can be easily 


justified. We then pass to the limit A — > +00 in order to construct and characterize solutions 


to (39). These points are addressed in the proof of Proposition 3.3 in Section pi and the 


proofs of Theorems 2.1 and 2.2 are then straightforward owing the previous propositions. 


4 Existence theory and estimates 


This section is devoted to the existence of solutions to (24), their regularity, and to the 


derivation of some important estimates that will be used throughout the proof. We introduce 
first the Green’s function 

_ /- : -- 

G^(z,x) = - —t— with x = (z,x) and = k^J 1 + ia £ e 4 /fc 2 . 


47t|x| 


We then recast (24) into the integral form 

PI - T efil = u Q>e , 

where 


T e u(z,x) = -^- - [ Gl,(z — z', x — x')v( — , — ^)u(z', x')dz'dx' 

47T 7(0,L)xK 2 V £ £ 7 


ua,e(u,z,x)= / L Qe (u,z,x,x')fa(u, — )dx' 
JR 2 v £ 7 


( 41 ) 
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and 


e iKV\ z ~ L s\ 2 +\ x -x '\ 2 / , n I 

Lo !£ (uj,z,x,x') = (z-L s ) — --— -—T [iK [\z-Ls \ 2 + \x - x'\ 2 ) 2 

Att(\z — Lg\ 2 + \x — x'\ 2 ) 2 V v J 

We have then the following lemma: 

Lemma 4.1. For all e > 0 and all k u E M, © admits a unique solution p ^ in 

W 2 ’ p ((- oo, L s ) x M 2 ) n W 2 ’ P ((L S , oo) x M 2 ) 



that satisfies (24) in the distribution sense and almost surely. 


Proof. First of all, since the potential V is bounded, it follows from Riesz compactness 
criterion (see [50 , Theorem XIII.66 pp. 248]) that the operator T e is compact in L P (M 3 ), 
1 < p < oo. Moreover, application of the Young inequality show that uq )E E L q (( 0, L) x M 2 ), 
1 < q < oo. The existence and uniqueness of a solution to (41) in L p ((0,L) x M 2 ) is then 
a consequence of the Fredholm alternative. Remarking further that uo e E L p ((—oo,Ls) x 
M 2 ) D L p ((Ls, oo) x M 2 ), the solution p £ u belongs to the latter space and satisfies (24) in the 
distribution sense. Standard elliptic regularity finally yields p £ u E W 2,p ((—oo, L$) x M 2 ) D 
W 2,p ((Ls, oo) x R 2 ). □ 


A first consequence of the latter lemma is that the Helmholtz equation is satisfied almost 
everywhere for z > L$ and z < Ls- A second consequence is that p ^ and d z p £ u both admit 
limits as z —> Lg. The Helmholtz equation then yields the following jump conditions across 
the plane z = Ls, 

P £ lA l S’ k ) ~PU L g,K) = /o(w,k) and d z pl(L^, k) - d z f u {Lg, k) = 0. (42) 


We then use these relations to solve the Helmholtz equation for z < Ls, for z E (Ls, 0), and 
for z > L. This will allow us to derive boundary conditions at z = 0 and z = L that will lead 
to some estimates on p £ u . We will need for this some properties of the A £ w defined in (28), 


which follows from some of the principal square root. For a complex number z = u + iv, with 
r / 0, the principal square root admits the expression 


\fz = — y= (\f-\/u 2 + v 2 + u + i sign(v)\J \/u 2 + u 2 — uj. 


(43) 


As a consequence, 

Re(A eia ,(«:)) >0, Vk E M 2 , and Im(A £iW (/c)) < Cyfafffi, for 1 — £ 2 |k| 2 /A: 2 > 0. (44) 
For the second inequality above, we used the fact that the square root is of Holder regularity 


1/2. We will also need the following expressions, that are consequences of (26)-(27): 


and 


ai(z,K) = 


b £ u (z,K) = 


y/\ e ,UK)e- ikuXc ’ u,{K)z/e 


X ^^ )z/£2 r 


PI(z,k) + - 


iku A £iU) (k) 


dzPl(z,K) 


- 7 


ik^j A E)tl j(/c) 


dzPl(z,K) 


(45) 


(46) 
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Solution for z E (—00 ,Lg) U (Lg,0). Since there are no sources at z = — 00 , there are 
no right-traveling waves in z < Lg and therefore 0 ^,( 2 ,«) = 0 for z < Lg. Moreover, since 
the medium is homogeneous in (— 00 , 0 ), the coefficients a^(z, k) and 6 ^( 2 , k) are constant in 


(—00 ,Lg) U (Lg, 0). Using the jump conditions (42), we find 


at(L+, K ) = Q ( UjK )e- ik »>**>W L s/*. 


(47) 


The coefficient is unknown at this point in (Lg, 0). We then eliminate it in the expression 
of Pu(z,k), which leads after direct manipulations to the following boundary condition at 
2 = 0: 


dzPU 0 ,k) + 


ik ul X £jU! (izj 


PU 0,k) = 


X £ ,cj(li) 


k), 


Vk E 


(48) 


Note that we used here the fact that d z p £ u (z,iz) and pf d (z,n ) are continuous in z, tz a.e. 


according to the regularity of Lemma 4.1 


Solution for 2 E (L,+oo). As in the previous case, there are no sources at 2 = + 00 , and 
therefore bf J (z,n) = 0 for z > L. Since a^(z,K) is constant in (L, + 00 ) and unknown, we 
can eliminate it in the same fashion as above to obtain the following boundary condition at 
2 = L: 


dzPKL, k) = p^L, k), V/v E M 2 . 


(49) 


We then use the boundary conditions (48)-(49) to arrive at the following result: 


Lemma 4.2. The wavefield p £ u satisfies the following estimates: 


OLfE 


112 


/ £e(A eiW (K))|p£(L,/«)| dn+ / Re(X £tU (iz))\pf jJ (0, /c)| dn-\ - ||Pwllz, 2 ((o,L)xi 

J r 2 Jr 2 u ' 

— C||/o( W ) Oil L 2 (R 2 )> 

f Im(\ £ ^(K))\p £ UJ (L,K)\ 2 dK+ [ Im(X E ^(iz))\pl(0,K)\ 2 dK+^\\VfJ\ 2 L 2 ((0L)xi 

Jr 2 Jr 2 u ’ ' 

< C\\fo(u), 0IIl2(r2) ^1 + > 

II PKL, OIIlrr 2 ) + ll^j(0) 0 IIl 2 (r 2 ) < Cjl/o^, -)IIl 2 (r 2 )- 

We have moreover ^ 

1™II^(V)||l 2 (r 2 ) < 2 ll/°( w ’ 0 IIl 2 (r 2 )- 


Proof. We start by multiplying (25) by p e u (z,K). Integrating in (z,k) over [0, L] x 


(50) 

(51) 

(52) 

(53) 
using 


boundary conditions (48)—(49), and taking first the imaginary part leads to 

.2 


CXrS 


Re(X £tU] (K))\p £ u (L,n)\ dn+ / Re(A e , w (#e))|p® (0,/c)| diz+—— ||^IIl 2 ((o,l)xR 2 ) 

JR 2 K'uj 

= 2Re ( / yjX Et0J (n)a £ UJ (L ^, n)p^(0, n)dn) . 


Since Re(A £ia ,(ft)) > 0 for all k according to (44), we then find 


^ Re(X E}U} (K))\p £ bJ (0,K)\ 2 dK 1 < 2 ^ \X £tU (K)\\a £ u (L^, K)\ 2 dK^j ^ (0, n)\ 2 diz 
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where K is the support in the k variable of /o- Since there are constants C\, C 2 and C 3 such 
that (the last inequality follows from (|44|) ), 

0 < Ci < Re(A £)W (/«)), |A e , w («)| <C 2 , e I^IMV«-("))We 2 < Vk € K, 

deduce from (|47|) that 


we 


Ik 


\p £ u(0,K)\ 2 dK<C |/o(u;,ft)| 2 fk, 


(54) 


and therefore (50). For the second estimate, we take now the real part, and obtain 

e 2 ' 

kui J(0,L)x 

_ ku) 

= 72 


\X7p £ U) (z,x)\ 2 dzdx + / lm(X £0J (K))\p £ u ,(L, n)\ 2 dn + / Im(A £ W (/c))|p* (0, n)\ 2 dn 
J R 2 J R 2 

f (l + e s u(-, -))\p £ u (z,x)\ 2 dzdx + 21 m ( [ , k)$j( 0, n)df 

J (0,L)xR 2 V v £ \J R 2 v 


Using the fact that Im(A £i£J (ft)) > 0, that V is bounded uniformly in (z,x), expression (47), 
estimates (54) and (50), the estimate (51) follows directly. Estimates (52) and (53) are 
obtained in a similar manner as (50), we instead multiply (25) by p^(z, /«)A ea ,(re) . In that 

case, we obtain 

IK(r,•)&(„*) + IIS< 0 ,•)&(„,) < ||/„^,.)e^l r ’"( J *.«™ i V= J || i2(R1) || # = (0> .)|| ia(R1) , 

and we conclude using Young’s inequality. This ends the proof. □ 


Following definitions (45)-(46), it is then direct to estimates a £ u and b £ u from Lemma 


4.2 


Corollary 4.1. The following estimates are satisfied: 


(l a e u (z, k )| 2 + | %{z, n)\ 2 ^Jdzdn < -%—\\ / 0 (w, OHla^a) 


'(0,L)x{|re|<|fc„|/e} 


ale 4 


'{M<|fcw|/e} 


1^(0, k )| 2 dn < C\\f 0 (u, •) llz, 2 ( R 2 )- 


(55) 

(56) 


Proof. The bound (55) is a consequence of (45)-(46)-(50)-(51) and the relation 
\/oiu,£ < |A £jW | < -(1 + |A £;(X ,1 2 ) < -(1 + a UjE + (1 + £ 2 \k\ 2 /k 2 )). 


For (56), we use (26) at z = 0, together with (50), along with the fact that y/2Re {\ £)CJ ("))> 
|A £jW (k)| when |/c| < k u /e, and the calculation below: 


/{M<|fc w |/ e } 


1^(0, n)\ 2 dzdn < C\ 


< C 2 


'{M<M/e} 

/ 

TM<MM 

2 


d«|A e , w (rc)| [\p £ u (0,k)\ 2 + |/o(w, k)| : 


dnRe(\ £}UJ (K))\plj(0, k)\ 2 + C 2 \\fo(u} 

> ') lli 2 (M 2 ) 


< C3 ||/o(^, Oil l 2 (m 2 ) ■ 


This ends the proof. 


□ 


Let us remark that the absorption term a £ is introduced precisely in order to obtain (55), 
which allows us to control a e u and b e u on (0, L ) and not just at z = 0 and z = L. The estimate 
is used in the proof of the fact that the coupling with evanescent modes in negligible, see 
section |HJ 
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5 Technical results on iterated integrals 

This section is devoted to crucial technical results that will be used throughout the paper. 
The following fact is of importance: for f(ui ,.... u n ) an integrable function, invariant with 
respect to any permutation a. that is 


f (^1 > • • * j Un) f (^G(l)) • ■ ■ > ^ a(n )) j 


we have 


/ f(ui,...,u n )dui...du n = — ( f(ui,...,un)dui...dun, 

JA n (z) n\ J[0,z] n 

where A n (z) is the simplex defined by 

A n (z) := {(iti,..., u n ) € [0, z] n , s.t. 0<Uj<Uj-i Vj € {2,... ,n}}. 


(57) 


(58) 


The next proposition shows that the random field V satisfies a long-range property in the z 
direction. 

Proposition 5.1. For all zq € M and (x,y) el 2 x M 2 , we have 


E[V(z + z 0 ,x)V(z 0 ,y)\ ~ ^ 

z—H-oo z 1 


Ro(x — y) with 0% := — ( [ uQ(u)e u ^ 2 du) 

27tW-oo ' 


The correlation function Ro is defined in (13). 

The proof of this proposition follows the lines of |25i Lemma 1]. We give its proof below 
as a preliminary to the proof of Proposition 5.2 further. 

Proof. Let us first note that 


E[V(z +zo,x)V(z 0 ,y)\ = f xn(dq)e y) E[Q{B^{z + z 0 , q))Q(Bsj(z 0 , <?))], 

Js 

so that we just need to investigate the term E[0(Sj-,(z + zq, q))@(Bsj(zo, g))]. The analysis is 
based on the Hermite polynomials defined by 

g( 0 (n) e - u2 / 2 

Hfiu) := (-1 )' 1 9 . , with g(u) := , (59) 

9{u) V2vr 

which form an orthogonal basis of L 2 (M, g(u)du): 

(Hi, H m ) L 2( KiS ( u ) du ) = t'dim- (60) 

Decomposing 0 with respect to this basis, we have 

0 ( U ) = E Wllere 0 * := ( H l’ Q )Ll(R,g(u)duy 


We will also use Mehler’s formula which, for two centered Gaussian random variables such 
that E[A 2 ] = E[X|] = 1, yields 


E [H l (X 1 )H m (X 2 )\ = l\E[XiX- 2 \ l 5 lm . 
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Then, we obtain 


E[Q(B^z + z 0 ,q))Q(B^z 0 ,q))\ = £ + z 0 , q))H m (B^z 0 , g))] 


i!m! 


Z,m> 1 

= E \r l ^{z)R l {q,q) 


i> 1 


= e;r s (z)+£f4(*>- 

l>2 L 

Moreover, following we have z^r l ^(z) —0 as z —> +oo for l > 2, and also 

£ |f rj,w| < c£ f < c<e,e> IS < +«>, 

l>2 l> 2 


for z large enough. As a result, using dominated convergence for series, we obtain 
z^{^Z + Z 0 ,q))Q{B^q))] z ~ +oq Cs) Sl 


which concludes the proof of Proposition 5.1 


□ 


Since the proof of Theorem 3.3 is based on a moment technique, we will be required to 
compute moments of the form 


1 

£-n(2— s) 


E 


’An(z) 


rt 

[n®( e 


3 = 1 


f Uj 

'hi £ iPj 


<Pe(z,Ui, . . . ,U n )d'Ui ... du n , 


where n is an even number (otherwise this moment is 0 by symmetry), and ip e is a bounded 
function. The following result is extensively used in the forthcoming sections up to simple 
modifications. It provides us with crucial uniform (in e) bounds as well as with an important 
convergence result. 

Proposition 5.2. For all even number n> 2 and s = 2 — Sj/2, there exists a constant C > 0 
such that 


sup sup ——r / 

£6(0,1) PI ,.:,Pn £ n ( 2 ~ S ) J[0 tZ ]r 


E 


[ n e(£h(->Pl))] \dui... du n < C n n n 


3 = 1 


and, for the Cf of Proposition \5A 

1 r 


lim 


s n(2-s) J An{z) 


E 


[ n @(^(— ,Pj)) <Pe{z,ui, • • .,u n )dui... du r 

3 =1 £ 


= lim CL 

£-5>0 P 


■n/2 


/ E II ^ Pa,P ^ tp £ (z,ui,. ..,u n )dui ...du r 

,A ^T {a L p)eT\ u “- u ^ 


where cp e is a uniformly bounded function in e. Here, the sum runs over the pairings J~ of 
{1 ,..., n}, and the limit e —>■ 0 is uniform with respect to (p\,... ,p n )- A pairing over vertices 
of {1,, n} is a partition of this set made of n/2 pairs of couples (a, /3), for which a < (3 
and such that all the elements of {1,... ,n} appear in only one of the pairs. Note that the 
number of pairings behaves like ?W 2 , which appears in the estimate above. 
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Proof. The proof follows some of the ideas of |M]. For the first result of the proposition, 
we decompose ©(A” 1 •) over the Hermite polynomials (with resulting coefficients @ n ,i), and 
obtain 


E 


[n 0 N 7 ^ 


= E[ll0(r^( 


3 = 1 


,Pj 


3 = 1 


0 


n,Z, 


E II ■ 

ij>i V-1 


'll 


3 = 1 


Note that we introduce the factor 


An :— 


(n — 2) 1 / 2 


in order to force the convergence of a series, as will be explicit further. This is a key point 
of the proof. We want to use now [531, Lemma 3.2], which states that for n > 2, and a 
(Xi ,..., X n ) mean zero Gaussian vector such that 


E[A 2 ] = 1 and | E[XjXi]\ < 1 V(j, l) € {1,...,n} 2 with j ± l, 


we have 


E 


n 


j=i 


hi ■■■In'- 

2 q (ql) 


E 


r *l jl r *2i2 ' ' ' r iqjq 


if l\ + ■ ■ ■ + l n = 2q and 0 < li,... ,l n < q 
0 otherwise 


(61) 


where rjj = E \XiXf\, and 

I Uli • • * i In) = {(*1> Jl) • • • j iqijq) £ {!)•■■) n } ^ i s -t. 7^ J/3 V/3 £ {1, . . . , (/} 

and all index r 6 {1,..., n} appears l r times}. 

Above such a Gaussian vector is said to be standard. Nevertheless, because of the factor 
An we cannot apply (61) directly, we first have to make use of the following multiplication 
theorem m : 


p/2] 


Hi(X n u) = Ajj ^(1 — A, 


— fc 


ll 


k =0 


2 fc (/ — 2/c)! A! 


Hl-2k(u). 


Specializing (61) to our case, we find 

.(«.(» 


e[ii^-= 

i=i 


ri!--u 

2<?g! 


E n 


= J(h,..j n )P=l 

if l\ + • • • + l n = 2q and 0 < h, ■. ■ , l n < q with lj := lj — 2 kj, 

0 otherwise. 

Let us remark that all the indices l are odd since 0 is assumed to be odd (0 n ,z = 0 for l 
even). Hence, lj = lj — 2 kj > 1 for all j = 1,..., n, so that q > n/2. Consider now the term 


Aq : n '■ — 


[0,zf 


f[ r * C m r Ujm )X(Pir»,PiJ 

m =1 


dui ■ ■ ■ du n . 
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We deduce from the definition of I(h,.. ■ ,l n ) that each of the ui,...,u n appear at least 
once in the product above. Keeping n/2 of them for integrating , and bounding r% by 
sup n |r^( , u)| = 1 for the others, and using the fact that r^(u) is even, we find 

/ \ q—n/2 ^ q / r z \n/2 

Aq,n < (2z) n/2 ^sup|r fl (©)|J ( sup \R{pi,P2)\) | r h(~)|d«] . (62) 

S. - ^ - .✓ 

<1 


We now need to estimate the cardinal of I(li ,..., l n ). For this, we use again (61) with 
X\ = ■ ■ ■ = X n = X where X ~ jV(0,1), and find, with now rj m j m = 1, together with (62), 




E[n^ 

3 = 1 




< c r 


du\... du r 


r z fu\ W 2 r JL 

i r*fc)r) E [|n^-»iW 


Moreover, we have 


E[|n^W|] <n(n-ip/ 2 ./rd. 

3= 1 


3=1 


according to [33J Lemma 3.1], which yields 


(63) 


E 


l[ 0 ,z]' 


[ n H h ( vs * 

3 = 1 




du\ ... du r 


< c r 


'0 


r Sj 


© 


du 


n/2 TT X " ljl 

\y 2 ]! 


\y 2 ] ! 


x Y, II (n-l) l ^ /2 - k ^X- 2 -lp - 

j=j =1 2 k ik j \J{l j -2k j )\ 

kj= 0,...,[lj/2\ 

After standard computations, we find for lj odd, 

y]{lj- 2kj)\ > 2fe/ 2 ©^(©/2] - kj)\, and (n - l)^/ 2 -^ < n 1 / 2 (n - l)fe/ 2 I- fc i, 
and then, with the binomial theorem, 

Vi/ 2 ] 


©/2]! 


^(n-l)© 2 ^(A- 2 -l) - 

2 A ' ; kj ! ^ {/, 2kX 

n l/2 Vi/ 2 ) 


< 


2 [lj/2\ 

n l/2 


v (n - 1 )Vi/ 2 ]-ki(\~ 2 _ I) k i __ 

] (n j ©!(©/2] ©)! 


£ ^(n + A n 2 -2)fe«. 


26 



































Hence, using again that all the indices lj are odds, we obtain 


l[0,z\ 


3 = 1 


< n n/2 C n 


du i .. .du 

(?) 




du 


n/2 n yj lev , I 

tt | vyn,Zj I 

i£i }b {h/2] [h/ 2]! 


(n + A- 2 -2)fc 


/ 2 ] 


< (A n ?r 1/2 ) n C" 
Now, let us consider 




(?) 


du 


n/2 


|©n,2Z+l| 

,Z>0 ^ 


E 

Z>0 


|@n,2Z+l| 


/[nM]-l +oo 

= E + E 

\ 1=0 l=[nM]/ 

:=I+II, 


10. 


n,2Z+l | 

7! 


where M is independent of n and will be specified later. In what follows, we just work with 
l > 1 since the bound is direct for l = 0 . For the first term, we perform an integration by 
parts in @ nj2 z+i using definition (59), and obtain 

®n,2i + i = A-^-l ) 2l f ©«(A ?u)gW(u) = A" 1 ©^, 


and according to (60), we have 

|©n,2Z+l| < \ 1 ||@ ( ' 1) ( w )llL 2 (R, 9 (u) ( Zti)l|-^ 2 z||L 2 ( ffiiff ( u ) du ) < X n 1 SUp |© (1 )(u)E( 2 0- 

u 

As a result, using that (2/)! < 2 2 \l\) 2 we obtain 

[nM] — 1 

I <C\+ C 2 n 1/2 Y, 2 Z < Cl + n 1/2 C nM . 

i=i 

For the second term II, we have after 21 integration by parts, 

©n,2Z+l = K 2 \~ l Y J @ { 2 °(A n l u)g {l \u)du, 

and therefore, according to ([9]) , using that l\ > e(l/e) 1 , 

|@n,2Z+l|_v ^ X~ 21 


E 


i\ 


l>[nM] ' l>[nM) 

Then, setting M > eC q, we have II < C. Hence, 


<C5: 0 (eC|)‘. 


l>[nM] 


'[0,zY 


E [ II 9000 ) 

3 = 1 


du\... du n < n n ! 2 C n 


rsj 


(?) 


du 


n/2 


We finally conclude by estimating the term involving r^: following ([8]), there exists z e such 
that for all z > z e , we have \r^{z)\ < C\z\~^, and therefore, for all z > z e , 


I* |du < C(e + e * jT \u\^du) < Ce 2 ^~ 2 \ 


(64) 
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with s = 2 — Sj/2. 

For the second result of the proposition, we decompose 0 itself over the Hermite polyno¬ 
mials to obtain 


n 

e [IWb*(t'»)) 


E 

l a > 1 

a^{l.n} 



n 

E [n^( s *( 7 ’»)) 


©1 E II r fi(~- —)R(j>°,Pp) + Rn( u 1 ,■■■,««), 

T (a,p)er £ 


with 


Rn( u 1. 


and 


= E E 

j = l Sj(ll,...,l n ) 


Sj(h, ■ ■ ■, In) = {h = 1 for k < j; lj G {2,..., n}, Ik 6 n} for /c > j}. 


According to (61), R e n can be recast as 


= E E E n 

i=l g>n/2+l S jiq (h,...,l n ) ™=1 


@lm 

u 


n 

E[n 


m= 1 


where Sj_ q (li ,..., l n ) = Sj(li ,..., l n ) fl {Zi + • • • + l n = 2q}. Let us emphasize the fact that 
q > n/2 + 1 since there is at least one index lj greater than 2 and n is even. This is what 
will allow us to gain some extra powers of e to obtain the convergence to the leading term. 
We need to estimate for this the term A q ^ n for q > n /2 + 1 in the same way as before. Since 
rsj and R are bounded by one, we directly find, for all ..., i q ,j q ) G I(h, ■.., l n ), 



n £ Ujrn )R(Pim,Pjm) \ du l ■■■dun < 


l[0,z\ 


n/ 2+1 

n 

m =1 






d'UL\ ... dun) 


where ■ ■ ■, ^ n / 2 +i>fn/ 2 +i) re P ea t J twice. Since n/2 + 1 is odd, only one other index, 

denoted by jf. appears twice. In that context, two cases are possible. In the first case, we 
have a term of the form r^iiuj — Uj')/e ) (if any there is only one), and 



< 


Cie 2 * if fle (0,1/2), 
C'[elog(l/e) if Sj = 1/2, 
C'/e if Sj G (1/2,1). 


Using then (64), 




— «.,/ 

Jm 


£ 


\du\...du n = (/ du J dv\rsj( -)| dudv^j 

x du dv ri ( - )dudv 

Jo Jo £ 

< io g (i/ e ). 


n/2—1 


If we are not in the first case, we have a term of the form r^{{uj — Uj')/s)r^((uj — Uk)/s), 
k ^ j'. Using then the Cauchy-Schwarz’ inequality with respect to Uj, a change of variable, 
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the fact that r$ is even, and again (64) leads to 


l[0,z] 


n/ 2+1 

n m- 


m Jm 


m=l 


)\dui ... du n < C ( I du J dv\rs)( - )|dudv 


n/ 2-1 


r 2z 


du ri - 


< C' e hn/2^ A ( 1 -i5)lo g (l/ e ). 

As a result, bounding the cardinal of I(h,.. ■, l n ) in the same way as before, we obtain 


/ sup 

'[0 ,z] n 


E 


m =1 


[II H u{B%(-yiPm)) du\... du n 

n 

Cs * n / 2 £ bA(i-fl) log(1/e) | E [ Y[ H lm (X) 


m =1 


and therefore, using (63), 
1 


/ sup 

E 


• 5 ^n) 

/[ 0 , z] n Pl, — ,Pn 





<^Mi 0 g(i/ £)E e e n 

i=l <j>n/2+l (/i,... ,; n ) m=l 


VU- 


i®iK»-1)' 


< £* a(, - s) log(l/e)nC" £ 

y>i vf- 

According to © , we have 

| 6 j| = J @ {l \u)g(u)d' 


< CC l 


so that the sum above is finite, and which shows that for n fixed, the error term R e n converges 
to zero as e —> 0. It remains to treat the leading term. For this, we write 


y_ 1 _ f 

j £ n(2 S) J^n(z) 


n w 

(«,/8)eF 


u n - 


up 


) 


e^csj 


\u a - UrV* 


< 


(n — 1 )!! r 1 


n 


! L e 2 ( 2 -s) 


'0 JO 


du i ... du n 




E^Cf, 


- n\f> 


u — V 


dudv 


1 n /2 


where (n — 1)!! = n!/(2 n / 2 (n/2)!) is the number of pairings of {1,... ,n}. According to ([ 8 ]), 
for any 77 > 0 and z e such that z > z e , we have \r%(z) — < rjcsj\z\~^, and as a result, 


1 


e 2 (s—2) J Q 


a: n (^) 




I u — Up 



u — u| ^dudv 

J |u- 

-v\>£Z e 

+ e I 

rsj(u — v)dudv 

J 

\u—v\<Z a 

+ 

/ | u — v\~^dudv 

J \u—v\<£Z a 


This finally yields, for all > 0, 
1 


,— 1 f z f z fu — v\ e^Cfi , , [ z f z , , p, , , 

lim-^7—grr / / rfi - — i- -T- dudv < tjcp, / / m — u dudv 

e^ 0e 2{ s -2 ) J q J q £ ) u-v* ~ 1 Jo Jo 


£—>0£^( S 2) 

The proof of the proposition is complete. 


□ 


29 












































6 


Proof of Proposition 


3.1 


Let 0 = (01)02) £ Cq°(M 2 ) x Cg°(M 2 ) be a test function such that supp(cj)) C B(0,r,p). In 
order to control the transition between propagative and evanescent modes, we introduce the 
following integer 

n £ :=mf(n>0 s.t. {Q n E Q” >e s.t. 1 < e|Q„| 2 /fc 2 } / 0) 
where Q'f e is the set 

Q^,s = {Qn = ~ Qi - Qn s.t. (fi,gi,...,9 n )e5(0,^)x5" £ («:)}. 

Above, 5 ” £ (k) is defined by 


S™,e( K ) '■= {(qi,---,Qn) e 5" s.t. |<2z| < IM/ e V/ e {i,...,n}|. 
Let us remark that with this definition, we have, for all £ < £o, 
l /\k \ \ 

- rj and \X £ ,oj{Qn)\ > rj := s/ 1 - e 0 , 

Ls v V e 7 


(65) 


( 66 ) 


for all Q n E Q” e with n < n £ . The integer n £ measures the number of iterations it takes for 
the momentum to be at least of order 1 /y/e. Note that this order is arbitrary, any order 
of the form e ~ a , a > 0 would work just fine. Integrating then (29) in z and iterating n £ — 1 
times this relation, we obtain 


dy( z , K ) 

&&(*>«) 


n £ -l 


= J2 du (n) m(dq (n) ) i2£(«j,Qj-i, Qj) 


a£,(0, Q n ) 

fe(0,<3n) 


+ / du (ns) / m(dq (n£) ) JJ H^(uj,Qj-i,Qj) 

JA ne (z) JSZ ? s (k) 


3 = 1 


Oy( u n E i Qn e ) 
bfj{Un E , Qn s ) 


:= I e (2;, At) + J e (z, At), 


where Hf, and A n (z) are defined by (30) and (58), m(dqhd) := m(dqi) ... m(dq n ), and 


Qj := k — qi — ■■■ — qj, and Qo := At, . Let us point out the important fact that iterations 
are stopped before evanescent modes appear in the series (that is before the second term 


in the r.h.s of (29) enters the expansion), and this gives rise to only the two terms I £ (z,k ) 
and J £ (z,k). The first one is the driving term, and the second one is a reminder. The 
latter will be shown to be negligible, the main argument being that it is essentially of order 
1 /n £ \ after integration over the simplex A Ue {z). Writing a complete series expansion, i.e. 
choosing n £ = +oo (which then sets J £ to zero), requires to handle the transition propagative- 
evanescent mod es which is more difficult than treating the reminder J £ . Note also that in 
we are only interested in the convergence of a £ u (z, k) and b £ u (z , n) for k in the 


3.1 


Proposition 

support of 0, and we can therefore only consider the above equation for |ac| < \k u |/e. Thus, 
we set I £ (z, k) = J £ (z,n) = 0 for |ac| > \k u \/e. We then have the following two lemmas, 
whose proofs are postponed to the end of the section. 

Lemma 6.1. For all z E [0, L\, we have, for all /i > 0, 

limp(|(j (z),4>) l 2 ( R 2 ) xL 2 ( R 2 )| > = 0- 


This first lemma shows that J £ gives a negligible contribution, and therefore that I £ is 
the leading term. In the second lemma below, we introduce an auxiliary process I £ that 
approximates I £ by letting the absorption a £ vanish in the complex exponentials and by 
replacing the X £ttJ in the denominator by one. 
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Lemma 6.2. For all z E [0, L\, we have, for all p > 0 

limP(|(/ £ (z) - / £ (^),</>} L 2 (k2)xL 2 (m2) | >H)=0, 


where 


n £ — 1 


I e (z,n) := 


du W 


— 0 J An(z) 


n=0 


/ m(dq (n) ) jj ff£(iij, Qj-i, Qj) 

'SS.e(K) j= 1 


af,(0, Qn) 

fc(o,Qn) 


where H £ is defined by (34), and I £ (z,n) = 0 for |re| > |fc w |/e. 
Now, with the notation 


L e (^, k) := 


we have 


(/ 

aUz) 


^w(*) 

V \ 

to. 


££(*)_ 


BUz,k) 


,<t> 


with L £ (z,n) = 0 for |k| > |/c w |/e, 


> h) 


L 2 (R 2 )xL 2 (R 2 ) 

< P(|(I £ ( 2 ) - I £ (z),(j)) L 2( R 2) xL 2( R 2)| > /V 3 ) 

(-^ ( Z ) — L (z),<t>) L 2( R 2) xi 2( IR 2)| > /V 3 ) + p(|(^ (a)i ^)l 2 (R 2 )xL 2 (K 2 )I > 


+ - 


and owing Lemmas 6.1 and |6.2[ it just remains to prove that, Vju ' > 0, 
lim P(|</ E (z) - ^(-)» i 2 (R 2 )xL 2 (R2) | > n') = 0. 


(67) 


This follows from the calculation below. After straightforward algebra and the change of 
variable k —»• k + q\ + • ■ ■ + q n , we have 


(/ (z) L (^)) ^ ) ) i 2 ( R 2 ) xL 2 ( ]R 2 ) 


n>n £ 


a £ (o) 

k(o) 


,Rt 


L 2 (R 2 )xL 2 (R 2 ) 


where R £ n {ifi) = 0 for |k'| > \k u \/e, and is given by, for |«'| < \k u \/e, 


KM ■■=(,, f. 

JA„(z) Js 


m(dq (n) ) Y[ H £ (uj,Qj-i, Qj) 
S£, e ( K ') j =i 


<M<2o) 

foiQo) 


! A n {z) 

where Qj = k' + q n + ■ ■ ■ + (?j+i, and 

«Sft, e («') := {(gi,..., 7 n) e 5" s.t. |Qj| < |M/ e , VZ G {1,... ,n}}. 

We need now to bound R £ n . For this, the following estimate, which is a consequence of the 
first result of Proposition 5.2 will be used several times in the course of the proof of the 
proposition: 


T — 

- L, n..F. • — 


l 


l n,e ■— 2n(2-s) SU P / 

£ ( > An) (n) J[0,L] 2n 


Qi q 2 


IL 

dn^\E[l[v(^,q 1 j)v(^i,q 2 j)}\ < (n 1 / 2 ^) 2 ™ (68) 

i=i 


Since Qo £ supp(cj)), the support of R £ n is included in K n = {|k'| < nrs + r^}, and then, 
according to 


E 


\R £ n\\L 2 (K n )xL 2 (K n ) 


< C' 2n E[|m(S)| 2n ] \K n \ (sup \fij{u)\) 2 T n , e /{n\) 2 

3 = b2 
uSR 2 

^ (n 1//2 C ) 2n . , 2 

< —^—{nrs + rfi) 2 . 


(n!) 


H2 
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Using finally (47) and estimate (56) in order to bound a^(0) and b £ u { 0), as well as the Markov’ 
and Cauchy-Schwarz’ inequalities, we find, 

H\(I ( Z )~L ( 2 )) ^)l 2 (R 2 )xL 2 (R 2 )I >A0 

^ ^Wfoiu, -)IIl 2 (r 2 ) y- / F r,|p £ ||2 n 1 / 2 

-■- zL l 111 W K n\\L2(K n )xL\K n ) 




n>n £ 


< 


c 


J2 -/=X nr S + r</>), 

n>n r V^ ! 


which concludes the proof of Proposition |3.1| We end this section with the proofs of Lemmas 

0 and El 


Proof of Lemma 6.1 The proof is very similar to the one above, and we only detail the 
differences. The main ingredient is the fact that J £ is proportional to l/n e ! after integration. 
After the change of variable n —>• k + q\ + ■ ■ ■ + q n , we have 


(J (z), ^) i2 ( R 2 ) xi 2 ( R 2 ) du r 


K( u n e ) 


i Hn e {Unf) 


L 2 (R 2 )xL 2 (R 2 ) 


with H n (u n , k') = 0 for |k'| > \k u \/e, and for |k'| < |fc w |/e, 

H n (u n , k ) := / du (n_1) / m(dq (n) ) H*(uj,Qj-uQj) 

JA n (z,u n ) 


<MQo) 

<MQo) 


Here, Qj is as before, and 

A n (>, ri n ) u n -i) G [0,z] n_1 , s.t. u J+ i < Uj Vj £ {1 ,... ,n— 1}}. 


As a result, using the Cauchy-Schwarz’ inequality, (68) with n = n £ — 1, as well as estimate 
(55), we find 


E 


n e ( s - 2 )c ne ~ 1 

K J ( z )’ ^)l 2 (r 2 )xL 2 (R 2 )I_ — ,J(n e - i7f E 


du 


/ o 


'{N<l fc UM 


ck(|a^,(u,K)| 2 + |6*(u,/e)| 2 ) 


( s - 2)-2 

< _ (jm 

aj*sj(n e - 1)! 


According to the bound from below for n £ given in (66), and the fact that n\ > Cn n+1 ^ 2 e n , 
we find 


c s-4 C* o n z 

_ _ < ^ e r (8+l/2-25-n £ /2)log(n £ ) 




a e 


This concludes the proof of the lemma. 


□ 


Proof of Lemma 6.2 The proof simply consists in sending the absorption a £ to zero and 
using the fact that £ 2 \Q n \ 2 /k 2 —> 0 as e —> 0 for n < n e . Let us for this introduce, for 
l € {0,... ,n}, 


if, n (z,K) := [ du™ [ m(dqW) 

M„(z) JS™,(k) 




l 


3= 1 J = i+1 


0^,(0, Qn) 
^,(0,Qn) 
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so that 


n E — 1 n— 1 

I £ (z , K) - I S {z, n)= Y.Y, (Jl,ni z i «) - ^l+l,ni z ’ 


K 


n= 1 1=0 

After the usual change of variable k —> n + q\ + • • • + q n , we have 


( h,n { z ) h+l,n( z )i L 2 (R 2 )xL 2 (R 2 ) 
with If n (K r ) = 0 for |k'| > \k u |/e, and for |k/| < \k u \/e, 


«L(0) 
bU 0 ) 


T e 


L 2 (R 2 )xL 2 (R 2 ) 




/ du (n) m(dqW) H^UjiQj-iiQj) 

I An (z) JS2,'W / = 1 

n 

x [-^5(^+1) Qu Qi+ 1) - -^S( u J+i) Qh Qi+i)\ II n‘(.uj,Qi-i<Qi) 

3=1 +2 


^i(Qo) 
*^2 (Qo) 


where Qj = k/ + q n + • • • + Now, we deduce from (43) and the fact that the square root 
is of Holder regularity 1/2, that for all |g| < \k u \/s, 

\M+,M) - < 5 (\/{i-^l9l 2 A2 ) 2 + < e - (i - 7«l 2 /*3)) < - 2 ^ c . (69) 


With (|69|), (|44|), the definition of n e and (|66|), we then find, for l < n e , 

K,u(Ql) ~ ^e,u(Ql ) 


1 


K,UQi) 


- 1 


< 


K,Ul(Ql)K^j(Ql) 


+ 


KM) ~ 1 


KM) 


< 


Ca, 


1/2 


£,i0 


r\ 2 


+ 


Cy/e 


■■= le, 


and 

e 


ik u (\ S ,u(Ql)-Xe,u(,Ql-l)ui/£ 2 _ pi(.K,u>(Ql)-K,w(Ql-l)) u l/(. 2k w) 


Jk„(\e,v(Ql)-\ r e , u (Ql)-(\eAQl-l)-*lUQl-l)))ui/£ 2 _ 1 


< C 


6 ife w (l?e(A £ , w (Q i ))-A^((5;)-(«e(A e , w ((5 i _i))-A^((5i-i)))« ! /£ 2 _ x 


+ 


-fc„/m(A e ,„(Q i )-^£,^( ( 3i-i)) u i/ e2 _ X 


< c«y 2 . 


The rest of the proof is now classical and follows from the same techniques as (671: we first 
apply (68) to obtain 

E 


T £ II 2 

L l,n\\L 2 (K n )xL 2 (K n ) 


< ( 7 2 + a e ) C' 2n E[|m|(S) 2n ] \K n \ ( sup |0»|) 2 T„, £ /((n!) 2 n 2{n ~ l ~ l) ) 

3 = 1,2 
ueR 2 

(7 e+0( £ )C 2n 2 

-(nr 5 + ?>) . 


< 


n\ 


and then use (47) and estimate (56) in order to bound a^(0) and 6^(0), as well as the Markov’ 
and Cauchy-Schwarz’ inequalities to arrive at 

F (K I ( Z ) ~ I (")> ^)l 2 (R 2 )xL 2 (R 2 )I > k) 


< 


leC 


n s — 1 n —1 

EEC 


E 


V n =1 «=0 
7e — C7" 


Ite 112 

\ L l,n\\L 2 (Kn)xL 2 (K n ) 


1/2 


< — Y] —F^ri(nr s + r>) 

» n=i vn! 


This concludes the proof. 


□ 
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7 Proof of Proposition 3.2 


The proof is based on iteration techniques. We only prove the second point of the proposition 
since it is the most interesting one and the first point follows from similar calculations. We 
start by writing as an infinite series obtained by iterating (36), that is 


V u%S z '^ = with 7 u’ £ (z,k) := 

n> 0 


(70) 


where, for n > 1, 


7 j’ n (z,«) := [ du("> [ m(dqW)H 5 ; 1 > 1 (uW,qW)^(Q n ), 

7a„(2) JSsJk) 


and 


H; 


1=1 


■u 


;, Qj— l, Qj )■ 




by (34), and 


H^i,! 


denotes the (1,1) entry of the 


Here, 5 " £ (k) is defined by (65), 
matrix H^. As we will see, the leading term in Tff ,£ is the one obtained by the product 
of the diagonal elements of the matrices H^(uj,Qj-i,Qj). Any other term involving an off 
diagonal component introduces an oscillatory integral leading to a vanishing limit as e —>• 0. 
This will be proved further. The leading term is therefore 


X^(z, k) := K’ n ( z , «)» X°’ £ (z, «) := 


(71) 


n> 0 


where, for n > 1, 




’ n (z,K) := f duW f m(dqW)( V( Uj /s, 

Ve2_sy J An{z) j S u e{K) VfJ, J 

(72) 


with 


GJ(u<">,q<">) := J 


1 = 1 

Before getting to the core of the proof, we present some technical results that show that the 
two series above are well-defined, and that expectation and limits can be taken term by term. 


Lemma 7.1. The series (70) is well-defined, and we have, for all <ft £ L 2 (M 2 ), 

+oo 


E 


+oo 


«:*(*>■ = E ( £ cc = E E [ccw.Tw 


n =0 


n =0 


and 


lim E 
£—>0 


+oo 


RJ) J = i: limE[< 7 ^-'(^), 0 ) i3(K=} ]. 

n=0 

Similar properties hold for & 

Proof. We only consider ©, and just need to show that 

1/2 


SU P ( : 

n>l eS(0,l) 


E 


\Vf^z)\ 


L 2 ( R 2 ) 


) 


< Too. 
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This follows from a direct adaptation of the estimate on R fj in the proof of (67) in Proposition 


3.1 which yields 


/~<2n 

E[||7^’ £ (z)|| i2(R2) ] < — (nr s + r fo ) 2 , 
where is such that suppfo C B(0,r^ o ). This concludes the proof. 


□ 


that 


Owing Lemma 7.1, it is not difficult to show that X^ is the unique solution to (38) such 

E 


ll‘^'a;(' Z ) II L 2 (R 2 ) — ^ll/o( W ) ')IIl 2 (M 2 )' 

We now proceed to the proof itself and write 

KkM) ~ K(M = 


n> 1 


where 

T ^’ e ( z , k ):= ( du(n) f m ( dq (n) ) 

)£C n jAn ^ Js S,eW m = 1 

with Iq = l n ■= 1, and 

Cn-={l (n) = (h,...Jn-i)€{l,2} n ~\ s.t. 3m 0 6 {2,..., n — 1} 

with l mo = 1 and l mo -i = 2}. 

The set C n is such that there is at least one contribution of the off-diagonal of H^. Following 
Lemma 7.1 we can study the series above term by term, and in particular, 


(73) 




n> 1 


£—^0 


2 )J 


11/2 


where, using the second point of Proposition 5.2 


limE 
£—^0 


litres lll 2(R2) 

( k u C]l\ 2n - 
V 2 


dn [ E[m(dqi" ) )m(dq^ ) )]0°(gi )ri )0O(Q2^ 


^ {\ K '\< nr s+ r f Q } JS n xS n 

X H Zl / -^,n( 1 i n) » 1 2 n) > < li n) >q2 n) ) n R{q a ,qp)- 

(il,l,-..,h,n-l)e£n ^2,n (a,/S)eJ r 2,n 

(h,l,---,h,n-l)GCn 

We will show that the limit is actually zero. Above, r is as in the proof of Lemma 
the second sum is taken over all the pairings T 2 , n of {1, 2} x {1,..., n}, and is defined by 
(37). Here, we have introduced the notations 

Wt”T“ , ,q(" ) ,q? ) ) = iim 


7.1 


and 




du[ 


( n ) 


/ A n (z) 


du { 2 n) 


n 

(a,/3)eJ r 2,r: 


£f a -(0 Qa-(0,1)> Qa)£fp_ (0 Qfi-(0,1)> Q p) 


I u n - UR 1-6 


35 


















with 


££ j,l{ l ,u,p,q) := 


giku) (-^e,c o 

(?)- a e , w 

(p)W e 

if 

(. 3 , 0 = 

(i,i) 

i iku (Ag ^ 

j(q)+A £l , 

o(.p))u/e 2 

if 

U, 0 = 

(1,2) 

_giku) (A £)l 

u;(q)+A e , 

u(p))u/e 2 

if 

(i, 0 = 

(2,1) 

e ~ik u (\ £ 

,oj(q) A £ 

: ,„(p))u/e 2 

if 

(j,0 = 

(2,2), 

Now, let us consider a 

pairing T 2 ,n 

and l ( j 


i n) e C n . Using 

the notation of (73), take also a couple (ao,Po) such that ao(l) = 1 and «o(2) = mo for 
instance. The other cases follow exactly the same lines and are omitted. Using the following 
relation 

gir(u—v) 


— u| ^ = Csj 


11-3 


-dr, 


where := T(2H — 1) sin(7rid)/-7r with H = (2 — S j)/2, we single out the pairing (ao, fio) and 
obtain 


T nW l(") Z 11 m 

1 T 2 ,nK l l , h iQl ) Q2 ) — c f> 1™ 


dr 


|i-fl 


du 


(n) 


du 


(n) 


A n (U 




n 


^U*_(0,1)da( CK (-^)’ '““j Qa-(0,1)) Qa)^_ f0 ,i 1 ,i l g(^(l)) u /3> Q/3-{ 0,1) 5 Qp) 


/3-(0,l).'/3' 


(a,/8)&F2,n 


«« - ^ 


X e 0 | ^U 0 -(o,l)iU 0 ( a o(l)) w «o, Qa 0 — ( 0 , 1 )) Qa 0 ) 
X e - ir ^o £f Po _ (01)>ho (A)( 1 ), «/) 0 , Q/ 3 0 —( 0 , 1 ), Qa> ) ■ 


Integrating by parts (with respect to the variable it ao ) the function e lUa o^ r+k ^^,^{Qa 0 )+\i:,u:(Q ao 
with antiderivative 

e m ao (r-+^(Ae,„(Qc 0 )+A e ,„(Q QO _(o,i)))/e 2 ) _ ^ 

i{r + k U} (X £}UJ (Q ao ) + X e!U1 (Q ao _(o,i))) / £ 2 ) 
we find, using dominated convergence, 


<C\ 


du^ n) [ du(, n) TT K-wal a 


-l( 2 ) 

(a,p)ejz 2 , n 



(a,P)^(ao,Po) 


r 

g“aO-(0,l)( r +L(A £ ,u((3o(|)+AE,u(QaQ_(0,l)))/ E2 ) _ J 

dr 

x lim / 

£- ^0 J 

i(r + k OJ {X £tUJ {Q ao ) + A £)a ,(<5 Qo _(o,i)))/ e2 ) 

Wi-S 




A n (z) 


duj/^ 


n 

(a,P)eT 2 , n 

(<*j/9)^(“o+(0,l),/3o) 


U n — Ul 


l-fl 


x lim 
£—^0 


s i« ao (r-+MA e ,UQ ao )+A E ,UQ QO _(o,i)))/e 2 ) _ 1 


i(r + k u) (X e ^ UJ (Q ao ) + A £;a; (Q Qo _( 0j i)))/e 2 


dr 


11— 


wit 


r := (ui,i,..., «i, mo -l, '“l.mo+l, • • ■, ui ,„), and where A n _i(z) is the same simplex 

(n) 


as 


(58) with now n — 1 elements. Let us remark that we are working here with fixed ' and 


so that 


A eiC j(Qa°) + A £)I j(Q ao _(o,i)) — 2. 


Therefore, together with the Markov inequality and dominated convergence, the following 


lemma concludes the proof of Proposition 3.2 


(o,i)))/ e2 ) 

5 
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Lemma 7.2. For all a ^ 0 and u / 0, we have 


r I e iu{r-a/e s ) _ -o 

lim / - -7-77 pi— r 

e-s>o J \r — a/£ s \\r\ L ~^ 


dr = 0. 


Proof. Let ^ > 0 and r/ > 0 be small parameters, and let us decompose the integral into 
three parts as follows: 


| e iu(r-a/e a ) _ -o 


dr = 


a 


+ 


^ | e iu(r-a/e s ) — 1\ ^ 


J | r a/e^llrl 1 ^ ' J\r—a/e s \>ii/£ s J ? 7 <|r— a/e a |</i/e s J \r— a/e s \<ri ' \f Cl/£ s ||r|^ ^ 

We treat the last integral first and make the change of variable r —> r + a/e s r to obtain 


r | e iu(r-a/e a ) _ 

I\r—a/e 3 \<rj \r - a/e s ||r| 1_fi 


dr = 


- II 


< \u\ 


|r |<?7 |r||r + a/e s \ l % 
dr 


dr 


l\r\<T) |r + a/^l 1 -^ 


s(l-i5) 


< m e' 


< 


dr 


l\r\<r) | M - £ S/ q\ l h 


For the second integral, we have, with the change of variable r —> r/e s , 


r I e iu(r-a/e a ) _ ]_| 

Irj<\r—a/£ s \<ii/ e s \f o/e s ||r|^ ^ 


dr = e s 


< 


r | e iu(£ s r-a)/e s _ 

I £ a T]<\£ s i'—a\<fi \£ s r aUrl 1 ^ 


dr 


rs(i-h) 


dr 


|a| - £ s ri) 1 ^ J £ a rj<\r—a\<fi \r - a| 
< Ce s(1 - fl) log(l/e). 


For the last term, we find, after the change of variable r —>• r/e s , 


| e iu(r-a/e s ) _ -o 


i \r—a/e s \>ii/£ s \f Q/e s ||r|^ 




dr = e s(1 - f3) 


< £ s(!—i5) 


| e iu(r-a)/e a _ 
1-iD 


'\r—a\>fL y CL\\T 

f dr 


dr 


f\r—a\>fL y Cl\\T 


1—' 


which concludes the proof of Lemma 7.2 


□ 


8 Proof of Proposition 3.3 


The proof is split into two steps. We start by constructing solutions to the fractional Ito- 


Schrodinger equation (39) in the Fourier form. We prove the announced pathwise regularity 
in W^(0,L) and in (0, L), and show that, up to a phase shift, the obtained solution 
is the Fourier transform of the solution in the sense of Definition 1.1. The second step is 
to prove the convergence of the process Xf. We will show for this the convergence of the 
moments of Xf and use a regularized process. 


8.1 Existence theory for the fractional Ito-Schrodinger equation (39) 


As explained in the outline, we construct solutions via a regularization procedure. The 
solutions will be written in terms of Duhamel expansions, since, as in the proof of Proposition 
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3.2 it will enough to check the term by term convergence to obtain the convergence of the 
whole series. We then introduce the process X A , defined formally as 


«) = <#!(«) + 5Z x A ’ n ( z ,«), 


n> 1 


where 


X A ’ n (z,n) = {ik u ) n [ du™ [ m(dq( n ))e*G4uW,qW)^0 ( g n) 

JA n (z) JS n 

-w(dr m , q m ), 


l(-A,A) n \r m \ H ~ 1/2 


(74) 


(75) 


with 


G„(uW,qW) = - — £ ( \Qm \ 2 - |g m -i| 2 )« m . 
2k, ' 


(76) 


ra=l 


Here, we use the same notation as in Section [7j Moreover, (w(dr,q)) q ^s is the family 
of complex Gaussian random measures, independent of m(dq), defined by (23), such that 
w*(du,qi) = w(—du,q), and with covariance function 

E[u;(chc, gi)u; (du,®)] =--i?(gi, Q 2 )d(u — v)dudv. 

7r 

We introduce as well the regularized standard fractional Brownian field B^ and its k— th 
derivative b^ k := d^Bfj (bjj := b^' 1 ), for q € S, 


Bh{u, q) := 


b^ k ( u ,q) : = 


IH(2H-1) f A e iru — 1 


C' 


£ 




zi/2 w (dr , g), 


/_A zr|r| 

IH(2H - 1) /- A (ir) k ~ 1 e iru , , , 




£ 


/-A r 


Note that b^' k is well-defined since 


E 


chi / |m|(dg)|6^’ fc (u,g)| 


< CE 


|m|(dg)^(g,g) 


du 


fA I I 2 fc 2^. 


1 20 


/-A r 


2tf 


zr < Ca- (77) 


We will use the notation B^°° = Bh, with Bu defined by (16). Let finally 

X- 4 (V>)(z, k) := ik^an 


m(dq) ( Z d ue -*(l«-9l 2 -M 2 W( 2 M^,( U) K _ q )dB%{u, q), 
Jo 


for A G [1, +oo], whenever it is well-defined, and let us remind the reader about the following 
notation 


A a{B^(q)) 


1 


sup \Dj_ a [B^\ t -(s,q)\. 


r(i - a) o < s <t<L 


When A < +oo, we wrote X A in terms of the fractional integral defined in (19). Since B 
has smooth trajectories, it follows from [36] Theorem 2.4] that the fractional integral is equal 
to the usual Lebesgue integral. 

Our goals are then to show that (X a )a>i forms a Cauchy sequence in the appropriate 
metric space, and to show that the limit satisfies (39). We start by addressing the path 
regularity of X A , and by characterizing it as a mild solution to a regularized fractional Ito- 
Schrodinger equation. 


H 
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Proposition 8.1. The series X A defined by (74) converges in C^(0,L), and we have 

Xj-(z, k) = 4>°(k) T ikuCrH [ m(dq) [ 2 ) u K 2k »)X„(u, k - q)dBu(u , q), (78) 

Js Jo 

with o‘fj = Cjj/(H(2H — 1)). Moreover, the trajectories of X A belong to C“(0,L) , and for 
all z > 0, 

1 , 


\X W (t)IIl 2 (r 2 ) — \\X U (0)IIl 2 (r 2 ) — xll/o^, 


')■ 


(79) 


The path regularity of X A will be useful to justify future calculations. The proposition is 
proved in Section 8.1.1 The next step is to recover the mild formulation of ([5]) in the Fourier 
domain by passing to the limit A —> +oo in ( |78| ). We need for this uniform estimates in A. 
They are consequences of next two technical lemmas that follow the ideas of [28]. The first 


one (Lemma 8.1) is based on the Garsia-Rademich-Rumsey inequality [IB] below: 

Garsia-Rademich-Rumsey inequality: Let p > 1 and a > 1/p. Then, there exists a 
constant Ca,p > 0 such that, for any continuous function / on [0,L], we have 


L r L 


\f(t)-f(s)f<c &tP \t-sr ~ 1 f f 

Jo Jo 

for all (t, s ) E [0, L], and with the convention 0/0 = 0. 


\f(u)-f(v)f 


o Jo \u — u|"P +1 


dudv, 


Owing the inequality above, the increments (in time) of q) can be estimated as 

follows: 

Lemma 8.1. Let 9 e (0 ,H) and A E [l,+oo]. There exist positive random variables 
(Zg (q)) q es, such that for all q E S, 

I Bfi(t,q) - Bfi(s,q)\ < Zf(q)\t - s\ H -° 

with probability one, and for all p > 1, 

sup supE[|Z/ 1 (g)| p ] < Too. 

A>iqes 

Therefore, the trajectories of Bjf belong to C^~ 9 (0, L) for all 9 E (0 ,H), and we also have 


sup sup E A? a (Bfj{q)) 
A> 1 q£S L 


< Too, 


for a E (1 — H, 1/2) , and all p > 1. 


The second lemma we need consists in adapting |28l Proposition 4.1] to our context, and 
provides us with estimates on T A {if). 

Lemma 8.2. Let A E [l,Too], k > 0, a E (1 — H, 1/2), and if E Wfi +2 (0,L). We have the 
following relations: 

1. For all 0 < s < t < L, 


\\ T Wit) Wi s )\\n k < K 1>Uik / \m\(dq)A a (Bfi(q)) 


x dr 


rs(l + r s )\\if{r )\\ Hk+2 T 


r )\\n k 


(r - s) c 


T 


dy- 


r) -if{y)\\n k 

(r — y) a+1 
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2. For all t E [0, L\ 


\\T A mt)\\u k + 


rt \\i A { m)-z A ms)\\^ 

(t - s) Q+1 

\—2a i ^-a\ 


ds < K 2 ,a,k J \m\(dq)A a (Bff(q)) 


x Jo dr ^ t ~ r ' ) za + r Q )[ll^( r )H«fc+2 + H^( r )H H k + J 0 d v " Y ^ r _y^l} nk 


3. The following estimate holds: 


\\l A m\i-*,C,H k < A\a, k / \m\(dq)A a (Bfi(q))[\\iJ>\\ atk + sup \\i/>(r) \\n h+2 \ ■ 

JS re[0,L] 


Here, (Kj,a,k)je{i, 2 , 3 } are non random positive constants. 

Note that there is a loss of regularity in the q variable since estimates on X A (/if) in Hk 
require if E Hk+ 2 - This is due to the the lack of regularizing effects of the Schrodinger 
semigroup: regularity in time has to be exchanged for some regularity in space. The proofs 
of the last two lemmas are offered below. 


Proof of Lemma \8.1\ Let 9 E (0, H) and p > 1 such that p/9 E N*. According to Lemma 8.1 
with a = H — 9/2 and p = 2/9 , we have for all 0 < s < t < L, 


- B H(s,q)\ < |i- s\ H e C 


? L f L | B A {u,q) - B A {v,q)\ 2 / e ^0/2 


/ 0 JO 


u 


— v | 2 ^ T 


a.s., 


where C is a deterministic constant. Moreover, 
E[\Z^{q)\ p ] < C p L {p ~ l)e 


■■=ZHq) 


L f L E[| B A (u,q) - B A (v,q)\ 2 P^] ,e /2 


10 Jo 


I It — v\ 2Hp / e 


thanks to the Jensen’s and Holder’s inequalities. Using that Bf/(u, q) — Bf/(y, q ) is a Gaussian 
random variable, we find 


E[| B A (u,q) - B A (v,q)\ 2p / 6 } < C p ^E[\B A (u,q) - B A (v,q)\ 2 f/ 6 

r A I p ir(u-v) _ i 12 p /e 

<C P:9 ' 


L J-A 


\ r \2H+l 


-dr 


R p/d (q,q ) 


< C P)d \u-v\ 2Hp / e . 

This gives the first point of the lemma. Now, since 

1 rl 




T(a) L (t — s) 1 
+ (1 — a) 


t \Bff(y,q) - B£(s,q)\ 


(y-s) 


2 -a 


dy 


we have, following the same lines as above, with now using a 9 E (0, H + a — 1) such that 
p/9 E N*, 


E 


A <C pAa E[\Z^gW}<C rA J<’' 2 (g. q ). 


which concludes the proof of the lemma. 


n 
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Proof of Lemma 8.2. According to the first item of |28l Proposition 4.1] and the fact that 


r |1 _ e *(M 2 -|K-g| 2 )(r-!/)| rL | \ _ e *(N 2 -| K -<?| 2 ) 2 //( 2 M | 


dy 


— 7/i Q +l 


{r - y) 


< 


dy- 


•Q+l 


- ' ' ' U 2k u J 0 y<* 

< C\q\\2n - q\, 


we have 


\\l A ^){t)-l A {^){s)\\ Hk <(JH \m\(dq)A a (B^(q)) / dr Cr s \\ |2 • -<#(r, ■ - q) \\n k 


T 


a - - g )bh 

(r — s) a 


+ a dy 


r,--q)-i’(y,--q)\\n k 


(r - y) a+1 


which gives the first point of the lemma after the change of variable k —> n + q and basic 
computations. The second point follows closely the proof of the first item of [HZS i Proposition 
4.1]. The last point is a direct consequence of the first one. In fact, taking s = 0 in the first 
relation yields 


\\ T W{t)\\n k < Ci,a,fc sup ||^(*)lk fc+2 + sup 
L te[o,L] te[o,L] 


and also 

rA 


n k + 


-if(s)\\ n k 


o (t - s)" +1 


ds 


\\i i '(m)-r A ms)h t 
< C , 2,o,t(i - s) 1-0 


sup 

tE [0,-Z^] 


H k+2 + sup 
te[o,L] 




\w 


(t ~ s) 


a+1 


ds 


This concludes the proof of the lemma. 


□ 


Now, using the previous two technical lemmas, we have the following proposition, proved 
in Section 18.1.21 

Proposition 8.2. Let 6 E (0,H — 1/2). The family (A/) a>i converges in probability as 
A —> +oo in C^~ e (0, L), to a limit denoted by X^, which is the unique pathwise solution to 


X h 


,(z, k) = </°( k ) T ikuOH f m(dq) f e 2 l K l 2 ) u /( 2fc ^)^’ aj ( u , K _ q )dB H (u, q). (80) 

Js Jo 


A corollary (proved in Section 8.1.3) of this convergence result is the following. 

Corollary 8.1. The process defined by ^!^(z,k) = e A K \^ z K 2k A k) converges in prob¬ 
ability as A —>• Too in C^~ e (0,L) to ^^(z, k) = z P 2kuJ >X w (z, k), which is the unique 

pathwise solution to 


^ui{z,k) = 4^(0, At) - 


t K 


2 k 


w Jo 


^cj(w, n)du + ik u an / m(dq) / ^(u, n-q)dB H (u,q). (81) 


It remains to address the last point of Proposition 3.3 and to show that we can Fourier 
transform to recover the fractional Schrodinger equation of Definition 1.1. It is just a 
matter of switching order of integration. The proof is given in section 8.1.3 for the sake of 
completeness. 

Proposition 8.3. The Fourier transform realize a one-to-one correspondence between the 
solution of (81) and the ones of (§• 
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Finally, let us also remark that according to propositions 8.1 and 8.2 
r}> 0 


we have, for all 


ze[o,L] 


Pf sup ||^(^)||z, 2 (K 2 ) - -||/ 0 (w, -)IIl2(r2) 


= lirn P( sup \\Xu(z)\\ L 2 (r2) - ^||/o(w,-)IIl2(R2) >t})=0, 

A^+oc \ z£[0,L\ y 2 ' 


> T) 
• A 


1, 


which yields the conservation relation for and at the same time for 4/^ and 


8.1.1 Proof of Proposition 8.1 


The first step consists in studying the regularity of each term Xff ,A in the series (74). This 
is straightforward: since fff E Uk, f° r all k E N, recasting then X A,n in terms of as 


X*' n (z, K ) = ( ik u ) n / duW / m(dqWyc.duWqW) JJ &£(n m ,<z m )0°(Q n ), (82) 

JA n (z) JS n 


m= 1 


it is a direct consequence of (77), the Jensen inequality and the Fubini theorem that X A,n E 


C 1 ([0, L],T-Lk). We also obtain the recursive formula below by permuting order of integration, 

X™’ A { Z , k) = ik u aH [ m{dq) [ due^ K ~ q ^~^ K ^ u ^ 2k ^ X^~ 1,A (u, k — q)dB^(u, q). (83) 

Js Jo 


IS Jo 

We show now that for all A > 1 


N 


._ \ ^ \>r, 


(84) 


n =0 


is a Cauchy sequence (with respect to N) in probability in C°([0, L], LLk) for all k > 1. We 
will use the following result (see [2 Theorem 3.9 pp. 104]). 

Theorem 8.1. Let ( E,d ) be a complete metric space. A sequence (X n ) n >o of E-valued 
random variable converges in probability if and only if 

Vr] > 0 and v > 0, 3no > 0 s.t. Vn, m > no, P (d(X n , X m ) > rj) < u, 

that is (X n ) n >o is a Cauchy sequence in probability. 

Let us first remark that it is enough to work on the event 

r L 

IS \J 0 

where M > 0 is arbitrary. Indeed, using the Markov and Cauchy-Schwarz inequalities, 


E = 


\ m \{dq)(j^ du\bx’ 2 {u,q)\ + |6^’ 1 (0,g)| j < M j , 


¥ Uo dU Is l m K^I^’ 2 ( W, ^l > M ) - Jf ( E J 0 du J s \ m \( d< l)\ b H 2 ( U ^)\ 


1/2 


< C , 

“ M 


according to (77), and a similar estimate holds for the term involving bff 1 ( 0, q) after an easy 
adaptation of (77). Hence, on the event E, we have \b^' 1 (u,q)\ < M, for all u E [0,L], It 


then follows from (75) and the latter bound, that 


K’ A \\c°w,W < {1 +n|S|)‘||/„(w,.)|| Hl 


(■ CM) n 
n\ 
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The division by ra! comes from the integration over A n (z) as in (57). This yields the desired 
Cauchy property of the partial sum (X^ n )n. Now, using 

ll^(^-<Jv)llcO([0,L],W fc ) < C [ du [ \ m \( d Q)\ b H(u,q)\\\X^ - X^ N \\ C 0([0,L\,H k ), 

JO J >s 


it suffices to pass to the limit N —> +oo in (83) to obtai n (f78j ). The fact that d^Xf? e 
C°([0, L\,% k ) for all p > 0 is a consequence of (|77|) and of (|78[). It remains the obtain the 
conservation of the L 2 norm. Since we just proved that Xff G (^([O, L], TL k ), we can write 


d z X^(z, k) = ikujdH / m{dq)e * (|k 9,2 M 2 W( 2 M&>l( 2j q)x£(z, k - q), (85) 

Js 


so that 


(d z X^(z),X^(z)) L2{R2) = ik u a H Jj dnm(dq)e * (|k 9,2 l«l a W( 2 *«)^ g ) 

x x£(z, « - «) € *M, 

since 5, m and w(dr, •) are assumed to be symmetric. As a result, we obtain 

d 


d II yi 

dz" X “ 


Miiiw = 2 Re ((dz x *' [z) ’ x “' < 2 >W>))=°- 


which concludes the proof of Proposition 8.1 


8.1.2 Proof of Proposition 8.2 


The first step of the proof starts with the following lemma, proved further in this section. 
Lemma 8.3. The family (X^)a>i is Cauchy in probability on the complete metric space 

0^oo(0> L),da j00 ). 

Owing the latter lemma, let us denote by X^ the limit of {X^)a>i in probability in 


(Wff(0, L), dn.na), and let us prove that this limit satisfies (80). Note that the last point of 


8.2 


implies that (X^)a >:i is also Cauchy in probability in C^ 9 (0,L). According to 


Lemma 

the second point of Lemma |8.2| we have 


and 


II 1 + °°(X U - X*)\\ aik < K ajW / \m\(dq)A a (BH{q))\\X 0J - X^\\ aM2 


\ 11 + 0 °(XA)- 1 A (X^\\ a , k < f \m\(dq)A a (B H (q) - B£(q))\\x£\\ atk+2 . 

Js 

For the first term, we have for all q and M > 0, 

p(||Z +00 (AL - x£)\\ a , k > n) < f(||AL - X *\\ aM2 > v /(MK atk )) 

+ F^J^m(dq)A Q (B H (q)) > m ), 


so that according to Lemmas 8.1 and 8.3 


lim P 


(||X +0O (AL 


A —^+oo 

For the second term, we find in the same way, 


— X, 


u )\\a 


,k > rij = 0. 


( 86 ) 


»(||X+°°(^) -l A (^)\\ a ,k >v)< P (||^|| Q>fc+ 2 > M) 

+ P (J \m\{dq)A a (B H (q) - B^(q)) > r//(A/AA,fc)) • 
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Here, for the first term on the right hand side, using that the convergence in probability 
implies the convergence in law, the mapping theorem [3 Theorem 2.7 pp. 21], and the 
Portmanteau Theorem |BJ Theorem 2.1 pp. 16], we have 

hmsupP(||^||a,fe+2 >M)< P(||^|| aifc+2 > Af), 

and 

lim P( fc +2 > M\ = 0, since P(||T^|| Q fe+2 <+oo) = 1. 

M—>+oo V / 

Second, we have 


lim p( [ \m\(dq)A a (B H (q) - B^(q)) > rj/(MK ak )) = 0, (87) 

A->+oo v Js ' 


following the proof of Lemma 8.5 below. As a result, combining ( 86 ) and (87) we obtain that 
X u is a solution of (80). The next lemma addresses the pathwise uniqueness of solutions to 


(80), and concludes the proof of Proposition 8.1 


Lemma 8.4. Equation (80) admits a unique pathwise solution in W^(0,L). 


The section is ended by the proofs of Lemmas |8.3| and 8.4 


l oj 1 ‘''■(jj'N > > 7 ?/3) 


( 88 ) 


Proof of Lemma [iOj Let rj > 0, v > 0, N > 0, and let us write 

P(& > 0 o(^ +B ,^) > v )< P {d a ,oo{x£ +B ,Xffi 

+ P(<oc(<^,<iv)> 7 ?/3) 

+ nd a , 0O (x£ N ,x£) > 77/3), 

where X^N 7S defined by (84). First, let M > 0 and k v such that J2k>k v 2 ~ k < 77 / 6 , so that 

krj 

\,oo( x£ n ,x£) > 71/ 3 ) < p(E \\x£ n - x£\\ a , k > 77/6 


k =0 


< 


+ 


kf 1 

(J2 \\ X w,N - x u\\a,k > v/6, [ M (dq)A a (B^(q)) < M 

k =0 JS 

P(/ \m\(dq)A a (Bfi(q)) >m). 
j s 


In order to treat the first term in the r.s.h, we introduce 

:= \K- A 




— + l 


/o (t - s ) 

and find, using inductively the second point of Lemma |8.2[ 

Jn,k(t) < K 2 ,a,k [ \m\(dq)A a (Bfi(q))\ f dr\\X^~ ltA (r) \\n, 

Js 1 Jo 

+ f dr((t - r)~ 2a + r- a )J n _ hk (r) 


<K- 2 , a ,k / \m\(dq)A a (BH(q))t 2a / dr(t - r) Za r J n -i,k+ 2 (r) 


tfc+2 


< 


K 2 , a ,k \m\{dq)A a (Bfi(q)) Pn(t)\\f 0 (u,-)\\ Hk+2ri , 
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where (fi n (t)) n is given by 


A,(t) := 1 and /3 n+1 (t) := t 2a [\t - r)~ 2a r~ 2a f5 n {r)dr. 

Jo 


According to [25] Lemma 7.6], we actually have 


Pn(t) = 


r(l 2a) r,l —2 a 


[t 1_2a r(l - 2a)] 


r((n + 1)(1 — 2a)) 

As a result, since /o is compactly supported, we have 

r * i n 

, „ K a,k,Ll s \'m\(dq)A a (B^(q)) 

sn,A 11 " / j ~ ' " " ' s " 


X uj' A \\a,k A r(l — 2a)||/o(w, OII'Hfc 


r((n + 1)(1 - 2a)) 


(89) 


Hence, using (89), we have for all A > 1 and N v sufficiently large (but not depending on A) 

k-q 

IP(Ell<^-^IU,A:>ry/6, [ \m\(dq)A a (Bfl(q)) < M^j = 0, 

i n J S 


since we have for the previous event 


lx r : 


|a,fc — 


(MK a , k y 


T((n + 1)(1 — 2a)) 


Moreover, according to Lemma 8.1 and the Markov’s inequality, it is clear that 

Co 


supp( [ \m\{dq)A a {B^(q)) > m) < -j, 
A>1 y Js t M 


so that finally 


supP (d at00 (x£ N ,Xf) > 77/3) < 
A>1 M 


In the same way, we obtain 


supsupP(J Q!OOi 0 (A'^,A^ + - B ) > 77/3) < 

A>1B>0 ’ v M 


For the remaining term of (88), we have first 

km 


> v/3) < P(E H*-® - X ' A 


k =0 


u,N v ^OJ,Nr,\\a,k > v/® 


k =0 


(90) 


( 91 ) 


< E P (ll<tf - > vMkr, + 1 ))). 


Second, according to (83), we have for all n G {0,..., N v — 1}, 

{x n+l,A+ B _ X n+hA ){Zj k) = jA+B, ( ^n,A + B _ X n,A )(Z] k) 

+ iku f m(dq) [~ due- i ^ K - q \ 2 -\ K \ 2 > l ^ ) X^ A {u^-q)d{Bj I +B - B$)(u,q), 

Js Jo 
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and using the second point of Lemma 8.2 together with (89), we find 

11 "1^77.-)-1,^4-|-5 1^77.-(-1,^4.11 

W^uj \\a,k 

* n 

< / \m\{dq)A a {B A+B {q) - B$(q)) ^ \\X^’ A \\^ k+a{n+1 _ m) 

JS m =0 

r C a d ~\ n—m 

X Ci t a,k / \m\(dq)A a (B^ +B (q)) 

1 ' Js 

r ” [c 3 ,a,k J s \m\(dq)A Q (B A (q)) 

< C W J s MMM V (,) - B„( q )) £ 1 - r((m + 1)(1 — 2a)) 


Ci, a ,k I \m\(dq)A a (B^ +B (q)) 


(92) 


Now, using the facts that 

\^-< N nWa, k >V/m v + l 


N„ 


<£ 


X. 


n,A-\-B 


- x r : 


x,k > +1) +1)) 


n=0 


as well as 


vn,A+.B _ -yn^A 

lla 


, fc > c 3 ) 


< 


n,A+B 


CJ I|Q!,| 


+ 


k>C 3 , Js \m\(dq)A a (Bff(q)) < M, J^\m\(dq)A a (B^ +B (q)) < M 
P( J s \m\(dq)A a (B A+B (q)) > M) + P( jf \m\(dq)A a (B A (q)) > m). 


we have according to Lemma |8.1| the Markov inequality, and 


X 


n,A-\-B 


- x 


■ u ii a,k > C 3 ) < E + p( jf \m\(dq)A a (B A+B (q ) - £$(<?)) > C 5) , 

m\(dq)A a (B A+B (q) - B A (q))}/C 5 , n 


C A ^ 
< — + E 
~ M 


L Js 


for all M > 0 and all C 3iU sufficiently small (independently of A and B). Gathering all the 
previous estimates, we find, combining (90) and (91) in (88) and taking M —y +oo, 

\m\(dq)A a (B A+B (q) - B$(q)) 


W(d ai00 (X A + B ,X A ) > r,) < - + C 5iN E 


L JS 


Then, the following lemma allows us to conclude the proof of Lemma 8.3 
Lemma 8.5. Let a G (1 — H, 1/2), and 


Bi +B (u,q) := I 


We have 


{A<\r\<A+B} ir\r 

lim sup supE A a (B y } +B (q)) 
■ 4 ^+oc B> 0 q£S L 


g iru _ J 

Zip' u>(dr , q). 


H- 


= 0. 


As a result, we finally obtain with 


lim supP {d a ,oo(X A+B , X A ) > rj) = 0, 

A-s>+oo B>0 


which ends the proof. 


□ 
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Proof of Lemma \8.5\ Following exactly the proof of Lemma 8.1 we have for 9 E (0, H +d—1) 
such that 1/0 E N*, 


E 


A & (Bi +B (q))\ < K~ e& L H +^-' e [l + 1 /(H + a - 1 - §)] 

rL f L E[| Bi +B (u,q) - B i +B ( v ^)\ 2/§ ] dudv f 2 


L J 0 Jo 




and since B^ +B (u,q) — B^ +B (v, q) is a Gaussian random variable, one has 


iA+b 


supE 

q&S 


A & (Bi +B (q))} < Kg - sup B}^(q, q) 
J ’ q&S 


rL r L 


L J 0 Jo 


a 


\ e ir(u-v) _ 1 |2 1/0 dudy -i 0/2 

-ar 


{A<|r|<A+B} 


I2H+1 


m 


- V |2#/0J ’ 


The proof is ended using the dominated convergence owing that 

/■ Ig ir(u—v) -||2 


/{|r|>u} |r| 2 ^ +1 


-dr < dn\u — v 


■2H 


□ 


We end the section with the proof of Lemma 8.4 


Proof of Lemma 8.4 Let us consider two solutions of (80) denoted by X u and X^, and let 


Jk(z) = \\XM ~ Xoj(z)\\n k + / ds 


I XUz) - x^s) - (XM) - X u (s))\\n k 

— el“ + l 


(z-s) 


Using the second point of Lemma 8.2 we find 


Jk{z) < Cz 2a [ \m\(dq)A a (B H (q)) [ dr{z - r) 2a r 2a J k+2 (r), 

Js Jo 

so that considering 


N 


J(z) = v Km ^J k {z), 

iV—>+oo — 


k =0 


which takes its values in M+ U {+ 00 }, we obtain by the monotone convergence theorem 

J{z) < Cz 2a ( \m\{dq)A a {B H {q )) /“ dr{z - r)- 2a r~ 2a J{r). 

Js Jo 

As a result, according to the particular version of the Gronwall lemma given in |28l Lemma 
7.6], we have J(z) = 0 for all z >0. □ 


8.1.3 Proof of Corollary |8.l| 


First, it is clear \I/ A converges in probability as A —> +00 in & (0 ,L) to 4U. Moreover 

according to (85), we obtain 


®u(z,k) = ^(0,k) - 


l \K 


2 r z 


2 k 


w JO 


K)du + ik w aH / m(dq) / ^(u, k - q)dB^(u, q) 


(93) 
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and we only need to address the convergence of the last term. Introducing 


K, A (i/;)(z,k) = ik u a H / m(dq) / ^ A (u, k - q)dB^(u, q), 

JS Jo 

for A e [1, +oo] where B^°° = Bh , we have 

\\IC A ^ A ) - JC+°°(tf w )|kfc < \\(JC A - JC+°°)(tf w )|| ajfc -f \\IC A ^ A - ® w )|| a>fc 

< | K\a H [ \m\(dq)A a (B^°°(q))\\^ Lil \\ atk 
J s 

+ \k u \c H j \m\(dq)A a (B^(q))\\4> A - 

so that proceeding as in the proof of Proposition |8.2[ we obtain that, for all rj > 0 

lim f(\\/C a (* a ) - /C +0 °(* w )|| a>fc > rj) = 0. 

A —>+oo \ / 


Therefore, we can pass to the limit A -A +oo in (93) and obtain that 4U satisfies the desired 


equation. Uniqueness follows from the one-to-one correspondence between (80) and (81) via 
and X,.. 


8.1.4 Proof of Proposition 8.3 


With the notations of the previous section, we just need to show that 

<£ + “(* w )(z),$ l2(r2) = (27r) 2 </C +00 (^)(z),^) L2(R2) 

where, for A 6 [1, -fee] and the notation W^°° = Wh, 

iC A (ip)(z,x) = ik^ [ ip(u, x)dWff(u, x), W^(z,x) = a H f m(dq)e~ iq ' x B^(u,q), 

Jo Js 

The proof consists in approximating Wh by W^ with A < -foo in order to have sufficient 
regularity to justify the calculations. We thus write 


<£ +0 °(^)(z),^) l2(r2) = <(£+°° - iC A )(^)(z),0) L2(R2) + <^(^)(^)» L2(R2) 


with 


E 


||£ +0O (tfa,)(z)-£ A (tf w )(* 


sup A a{W H (x) - Wh(x )) 


< E 


< C\\fo{u, -)IIl 2 (r 2 ) supE A a (B%°°(q)) 
ges 1 


which converges to 0 as A -A- -foo according to Lemma 8.5 Following (77) and the fact that 
€ C°([0, L], H k (R 2 )) for all k G N, we can write, 

iC A {^>){z : x) = iku&H [ m(dq) I du^ u] (u, x)e~ iq ' x b^ I {rL,q), 

Js Jo 

and thus use the Fubini and Fourier-Plancherel theorems to arrive at 

<£ A (4' aJ )(z), < /.) i2(R2) = (2vr ) 2 </C a (^)(^), 0) l2(r2} . 

The proof is ended by passing to the limit in A by adapting what was done above for IC A . 
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8.2 Convergence of Xff 


We investigate in this section the limit in law of (Xf, (z),... 




The tightness of this 


family is addressed in Proposition 8.4 below and is the straightforward consequence of the 
conservation of the L 2 norm of Xf. The characterization of the limit in distribution of the 
family requires more work. We will use for this a moment method and the regularized process 
X£ for which formal calculations are justified. In order to fix the ideas, we will investigate 
first the moment of order one in Proposition 8.5, and then generalize to moments of any 


order in Proposition 8.6 
the section. 


The various results are put together in a conclusion at the end of 


8.2.1 Tightness 


This section is devoted to the tightness of the family (X£ (z),... ,X£ (z)) e . We have the 
following result: 

Proposition 8.4. For all z G [0,L], the family (X^z ),..., X£, M (z)) e is tight in L 2 (M 2 ) 
equipped with the weak topology. 

Proof. It suffices to show that the family of complex-valued random variable {{X^ (z ), 
is tight on C M for all cj)j 6 L 2 (M 2 ) (j G {1 ,... ,M}), which amounts to prove that 

M 

V 77 > 0, 3pi > 0 such that limP (EK^M.^I 2 >/•)<»■ 

3= 1 

This is a direct consequence of the following lemma: 


Lemma 8.6. We have, for all z G [0,L], 

ll^w(«)IU a (R a ) = II<?Ol 2 (R 2 )- 

The proof of the lemma is left to the reader. Since Xf has sufficient regularity (i.e. at 
least C 1 ([0, L], L 2 (M 2 )) almost surely), it suffices to adapt the proof of (79) to obtain the 
result. □ 


8.2.2 Moment of order one 

Proposition 8.5. For all z G [0,L] and <j> G Cg°(M 2 ), we have 

linrE[<^(z)» L2(R2) ] = ^lim^E^^),^^^. 

Proof. According to Lemma |7.1[ it suffices to show the term by term convergence of the 


series defining Xf. Moreover, since the integrand in (72) is L 1 in all variables, we can invoke 


Fubini Theorem to permute order of integration. Using the second point of Proposition 5.2 
we have for n = ‘In! (if n is odd the limit is 0 ), 


Iim o E[(^> 2 " , (z),^)] = (zfc w ) 2 "' 


dK(t>{n) 


du( 2n, )E[m(dq( 2n '))] 


x e 


J J ^2n'( Z ) JSn 

iG^uWqWUO { Qn )C$Y, II 

^2 n' / 


\u a - up p 


where the sum runs over the pairings F-in' of {1,... , 2 n'} and G n is defined by (76). We 
want to relate now the term above with X^ ,2n . We use first for this the Gaussianity of the 
measures w to find 

2n ' rA „irmUm , 


ge n 

T 2 n' 


R(q a ,qp) 


. _ = lim E 

\Uot — 14^31 A —^+oo 


[n 

m= 1 


'—A \r r 


I H 


—^w{dr m ,q m ) . 
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Moreover, since on the one hand, 


e[ii 


^ '‘A 


m=l 


/-A r„ 


iff—1/2 ^^Tmi Qrn) 


< C\u\ — U 2 


2H—2 


He (1/2,1), (94) 


and on the other hand that the integrand in (82) is L 1 in all variables, we can invoke both 
the Fubini Theorem and dominated convergence to obtain, for all n > 1 

limE [{K' n (z)A)} = (*M" [ d^W) f du™ [ E[m(dqW)] e ^»(« ( "W">)0O (Qn) 

J JA„(z) Js n 


£—>0 


X 


lim E [ IT 

—>+oo L 


n -A e irjuj 


A—>•+ oo 

= lim E [(Xu’ n (z),<j>)]. 

A— H-oo 


3 = 1 


I—A I r^\ H ~ 1 / 2 


—rj^w(drj, q-j) 


Above, (jff = lim £ is given by (37). It just remains to show that limit and expectation can 
be taken term by term in the series defining X^f. This is the object of the next lemma. 

Lemma 8.7. We have for all z £ [0, L\ 

Urn E [(X^(z),<f>)] = Y I™ E K*«’"(*). </>)] < 

—— A^-+oo 


A—>-+oo 


n>0 


Proof. This result is just a consequence of the fact that 

Y SU P ^[\\ X w’ n ( z ) lli 2 (R 2 )] 1/2 < +°°- 


n> 1 j4 ^ 1 


As in Lemma 7.1 

sA,n( ||2 


we 


have using 


E 




< k„ n / d« 


/S n xS n 


E[|m|(dq{ n) )|m|(dq^ n) )] 


' A n (^) 




(n) 


'A n (^) 


du^' 0 


r ” f A e ir ™ u ™ i i f A e ir ™ u ™ 

E U1 / , | r l \H-l/2 W ^ dr m^m) , 2 |g _ 1/2 w ( rfr LgT 

m —1 J —A \! m\ J — A | ' 777 , | 


A e * r m“m 


I /o(w, Ql,n)/o(^? Q2,n) \ 

_ c 2 

fo } ( n !)2 

2 C n 


< (nrs + r t ) 2 Yv? ( 2n _ 1)" 


Ld(0, 2 )2 


lui — W2| 2// 2 du\du2 


<(nr s + r fo ) —. 


Above, we used (94), the term (2n — 1)!! = (2n)!/(2 n re!) is the number of pairings of 
{!,... ,2n}, and the term (n!) 2 is a consequence of (57). □ 


This concludes the proof of Proposition 8.5 


□ 


8.2.3 Arbitrary Order Moments 

In the forthcoming computations, all indices with the subscript 2 correspond to the complex 
conjugate terms. 


50 


























Proposition 8.6. We have for all z E [0,L], frequencies (wgi, ..., , w 2 ,i, • • •, ^ 2 ,m 2 )j 

and test functions ■ • •, 4 * 2 , 1 , ■ ■ ■, 4>2 ,m 2 ) * n Cq°(M 2 ), 


Mi 


m 2 


j i =1 


liniE^ (^i^ 2 ) 7 <4,31 ) L 2( R 2) n (^2j 2 ( 2: )^2,j2) L 2( R 2)_ 

32=1 

Ml M 2 _ 

H+oo^t n ^ ui <n ^ >i >i i )l 2 (k 2 ) n ,J 2 ( Z )) ^2,32) 1,2 ( R 2) 

P2 = l 


A—>-+oo 

Proof. Using , we have 


31=1 


< +oo. 


Mi 


Mi 


M\ M2 +00 +00 


II (^ijj (")> ^1,3'i)l2(K2) H (^ 2 j 2 ( z )> ^2 j 2 )l2( R 2) E E E E X n : E X n’ 

31=1 3 2 =1 31=132 = 1 ni.j’i =0 n 2j2 =0 J n 


where 


,-ni-n 2 

X e = — _ 

n £ n(s- 1/2) 


Mi M 2 

n n dK ^32 k ^l7 2 02,32 O 2 . 


32 J 


31=1 


Mi 

n 

31 
M 2 

n 


■= 1 ^" M l (KlJi) 


m (^q£i w) ) 


32 

Mi 




m(dq 2 .; 


( n 2,3 2 L 


32 


^"ij'i ( z ) 


^ n 2.in ( 2 ) 


32 = 1 


du (n2 ’ , '2 ) 

aU 2,32 


.„ r , ("l.i, ) (»1.7i), “ 2 (™2.1o) ("2. 10 ). 

*ru ,G ” (uu > ^ n ^ 

31=1 32=1 

Mi n i,3i M 2 n 2,j 2 _ 

><n n ^(lll,3l,TOiji /£) 91,31, mi jj ) n n V(u 2 

,32, m 2,3 2 / £ 7 92,3 2 ,m 2 j 2 ) ) 


31=11711,^=1 


32 = 1 7712j 2 =1 


with 


Ml M 2 

ni := E niji, n 2 := E n 2j2 , and n, := n\ + n 2 . 

31=1 32=1 

As before, we need to show that limit and expectation can be taken term by term. 
Lemma 8.8. The series * s well-defined, and we have 

Mi M 2 


E 


[ II ( , ^1,3 1 ^)’^ 1 -3i)l2( R 2) H ,j 2 ( z ) >02 ,3' 2 )z, 2( R 2)] = E E t X n] 


as well as 


31=1 


Mi 


32 = 1 


M 2 


& E [ _II ( ;k ’-l,3i( Z )^bil)L2 (R 2 ) II (*W2, ia (-) >2,32 > L 2 (r2) ] =Ej^ E [ X n]» 


31=1 32 = 1 

Proof. As usual, it suffices to show that 


E SU P ^[l^nl 2 ] 1 ^ 2 < +°°- 

J n ee(0,l) 


(95) 
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Adapting once more Lemma 7.1 and the first point of Proposition [53] we have for s = 2—S)/2 , 


Mi 


E[|X n | ] < j ■ ■ ■ J 0?Kij 1 l{| reiJi | <nij . ir . s _|_ r .^} J^[ dK2j 2 l{|K2j 2 |<n 2 j 2 r s +r^ n } 


J 1 = 1 

|| C ni ’h || C n ^2 

iii n hw ! iii n 2 J2 ! 


m 2 

n 

J2 = l 


which gives (95). 


□ 


The latter lemma can be directly generalized to moments of X A using the fact that the 
measures are Gaussian. Using now the second point of Proposition [53] we obtain for n = 2n', 


limE[X*l = i ni ~ n2 
£—>0 


r. Ml M2 

■ / n n 

J J 31 =1 J2=l 

Ml M 2 r. 

**[n il L lu . 

L jl=li2=1 ^5 1 ji (ki,u)x5 2j2 («2j 2 ) 

/• / , ("1,1!) ("1,11 ) s 
X / dul^e 10 ^! ' q i.n ) ^(Qi Jl ,n 1 , 1 ) 

f , ( n 2 j 2 ) —iG,i(u 2 2 ’ 22 ^.q 0 2 ^ 2) )^077S y 

X / du 2,h e 2 2 V£(Q2ja I n 2 , J - a ) 

An 2,7, ( z ) 


n 


R{q a ,qf}) 

\u a - Up^' 


J~n (o£,/3){zJ ~n 

where the sum runs over the pairings J-" n of 

In := | (b Jij m i,ji ) £ {1)2} X {1,..., ALj} X {1, . . • , Jljjj} | • 


Moreover, in the same way as in Proposition 8.5 

R(q a ,q q0 ) 


ge n 


J ~n ((X,/3)E.J~n 


I u a - up ^ 


= lim E^ jQ [j 

j !=1 miji 


Ml n l,Jl ~j\ 'i‘ r l,ji,m 1 j ’Hl,ji,m 1 . 


A- 5-+00 L " . J-A In 1 . 

=1 mi j, =1 1 J-iJli" 1 !,. 


n 1 


II/—1/2 


~ w (drij 1:midl , f/lji.nujj ) 


M 2 n 2.j 2 e ir2 1 j 2 ,m 2 j 2 u 2,j 2 ,m 2 j 2 

X II II / , | r - IH-1/2 ^( rfr 2, j2 ,m 2 , J2 )g2J 2 ,m 2 , J2 ) 

j 2 =lm 2 ,i 2 =l A I "2 ,j 2 ,jti 2i j 2 I 
so that the proof is concluded by dominated convergence and the Fubini Theorem. 


□ 


8.2.4 Conclusion 

We have now everything needed to conclude the proof of convergence of X^. Consider first 


the limiting process X u solution to (39). Thanks to (40), the moment generating function of 
the random variable Y(z) = J2 jL 1 ■ ( z ) > 4>j) is perfectly defined for z fixed in [0, L], so that 

the law of Y(z) is uniquely defined by its moments. Then, since X A converges in probability 


to Xu according to Proposition 8.2, and since all moments of Y A (z ) := (Y A .(z), (j)j) 


(test functions 4>j in Cq°(M 2 ) are sufficient by density) converge according to Proposition 8.6 
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8.6 


m 


they necessarily converge to those of Y(z). Furthermore, since the limits of the mom ents 
of Y £ (z ) := J2jLi(Y £ .(z), (ftj) are the same as those of Y A (z) according to Proposition 
we conclude that the moments of Y e (z) converge to the moments of Y(z). Proposition 
finally implies that (X £ t (z),..., Xff (z)) e converges in law in L 2 (M 2 ) equipped with the weak 
topology to (X Ul (z), ■■■■, *u M (z)). 

Finally, convergence in law in L 2 (M 2 ) for the strong topology is obtained thanks to Lemma 
8.6, the Skorohod’s representation theorem [§1 Theorem 6.7 pp.70], and the following relation 


J™ KWW&pP) = Jm 


£— 


= jll/oCw, 


2) — ||'^(z) 11 £,2(^2). 


This concludes the proof of Proposition 3.3 


9 


Proof of Theorem 


2.1 


The proof is a direct application of Propositions |3.1| and |3.2| and Theorem |2.2[ As already 
mentioned in Section [3j owing the convergence results of Propositions |3.l| and |3.2| it is enough 
to check the convergence in law of 

■= JJ e- lut e lK ' x ^l(K)X^(L,K)dudK. 


Note that we cannot directly use here any continuity arguments of the map (X £ ) u e-)- (p| L )t.,x 
for the appropriate topology since we only previously obtained pointwise information about 
X £ in the variable u>. 

The proof is then done in three steps. First, we prove the tightness of l the space 
C°([—T,T],L^(M 2 )) for all T > 0, where L^(M 2 ) stands for the space L 2 (M 2 ) equipped with 
the weak topology. Second, we characterize all the accumulation points using the convergence 
of the moments of X £ . Finally, we obtain the convergence in L 2 ((— 00 ,+ 00 ) x M 2 ) and then 
in C°([—T,T],L 2 (M 2 )), where L 2 (M 2 ) is equipped with the strong topology. 

We will use the notation 


■■= (p|,l(V)» = (2tt) 2 (pi )L (t, •),<£), 

with 

= J e- wt ifi,(K)XZ(L,K)du. 

According to US Theorem 7.3 pp. 70], tightness is a consequence of the following lemma: 
Lemma 9.1. We have, for all 0 E L 2 (M 2 ), 

lim I5nP( sup \p e 2L (t, (j>)\ > M) = 0, 

Ad —£—>0 ' ££ rp 'y'j ' 


and for all rj > 0, 


lim limP( sup |p| L (*i, <j>) - j% L (t 2 , </>)]> rj)=0. 

T—>0 £—>0 \ tl -t 2 \<T 


Proof. This lemma is a direct consequence of Lemma 8.6 Let indeed <fi E L 2 (M 2 ). For the 
first point, we have 


limE[ sup |p2 iL (t,0)| < C'II0IIl2(r2) / du limE[|| || L 2 (r2) ] < + 00 , 
te[-T,T\ J J £ ^° 
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and the conclusion follows from the Markov inequality. In the same way, we have for the 
second point, for all r > 0, 


lirnE 

£—>-0 


SU P -P2,l(*2,' 

\u 1 -U2\<T 


<Ct I dw|w|limE [||^(L)|| L 2 (r2) 


which concludes the proof of Lemma 9.1 


□ 


In order to identify the accumulation points, we consider the finite-dimensional distribu¬ 


tions of p| L . We remark first that thanks to Lemma 8.6 


sup Hp 1 ,lWIIl 2 (r 2 ) 

te[-T,T] 

is uniformly bounded in £ by a deterministic constant, and therefore that the finite-dimensional 


distributions are uniquely characterized by their moments. Using once again Lemma 8.6 


in order to justify the use of the Fubini Theorem and dominated convergence, it follows 
from Proposition T3 that, for all M € N*, e [ -T,T ] M , G 

(L 2 (M 2 )) m , we have 


M 

iimE[ pl AtmAm) 

£—>•0 L , ’ 

m— 1 


= E 


M 

[ II PL^r 


0 


m= 1 


In order to go back to the original pulse pf defined by (20), we remark first that (35) holds 
for all cf) 6 L 2 (M 2 ) thanks to the bound (52) and the density of Cq°(M 2 ) in L 2 (M 2 ). Hence, 
using Proposition 3.2, we obtain the convergence in law of p e L to p® iu C°([—T, T], L^(M 2 )). 
To conclude, we use the Skorohod’s representation theorem [HI Theorem 6.7 pp.70]: there 
exist a probability space (f2,T,P) and random variables p s L and p^, with the same laws as 
p £ L and p° , respectively, and such that 


lirn sup |pi(i» -p° L (t,<f>)\ = 0 
£ ^°te[-T,T] 


— a.s, 


for all 4> £ L 2 (M 2 ). A direct consequence is that p^ converges P-a.s. to p° L in L^((—oo, +oo) x 


, since using 


one has 


1, 


1 ™IIPlIIl 2 ((oo,-Hx))xR2) < -||/o||l2(( —oo,+oo) X 


e-XT’ i-j mj-j u—» - — J — / - 2 1 

and the unit ball of L 2 ((—oo, +oo) x M 2 ) is weakly compact. Moreover, this convergence also 
holds in L 2 ((—oo, +oo) x M 2 ) with the strong topology because of the conservation relation 

2 II/o||l 2 ((-oo,+oo)xR 2 ) IIPlIIl 2 ((-oo,+oo)xR 2 ) |[PlIIl 2 ((oo,+oo)xI 

As a result, using the Plancherel theorem 

I e := / dudx\p £ L (uj,x) - ^> u (L,x )| 2 


£— 


0 


— a.s 


where 


Pl(L,x) = ±J e^pl(t,x)dt 


and 


&u(L,x) = ^ J e iUjt p° L (t,x)dt. 


Since /o(w,«) has a compact support with respect to u>, so do prj and according to (52) 

<CL. 


and (22). The Jensen’s inequality then yields 


sup ||Pl(V) -P°(t, 

te[-T,T\ 

This proves the convergence in C°([—T, T], L 2 (M 2 )) and concludes the proof of Theorem 
since almost sure convergence implies convergence in law. 


2.1 
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10 Proof of estimate (21) 


We use here the notation of sections [l] and [3j The core of the proof is the following lemma: 
Lemma 10.1. We have the estimate, 


C 


IIl 2 (R 3 ) < —=\\fo{u, -)IIl 2 (R 2 )- 


Proof. We already have an estimate on (0, L ) according to Lemma 4.2 which is 

C 


dlL 2 ((0,L)xR 2 ) < 


S\/OLi 


: ll/o(w,')|| l 2 (r 2 )> 


so that it remains to treat the domain z ^ (0,L). Consider first the case z > L. Plugging 
z = L into (45) and (46), and using the fact that b^(L, At) = 0, we find that p £ u reads 

PI(Z, At) = e ^XsA^~L)/P pULj K)> z>L 

We need to split the domain of integration in At in order to obtain appropriate estimates. 
Suppose first that £ 2 |At| 2 /fc 2 > 1 (evanescent modes), then Im(A eA ,(At)) > C^/a EtU} by (43). 
Together with (52), this yields 


r f 

dz / dn\p £ u (z,n)\ 2 < -=||/ 0 (w, 

4{£ 2 | K | 2 /fc 2 >l} s/OL e 




When 6 2 |k| 2 /A; 2 < 1 (propagative modes), consider the strictly decreasing function fi{x) = 
y/x 2 + b 2 — x, for x > 0. It satisfies fb(x) > fb{ 1) > Cti 2 for x G [0,1], 0 < b <C 1 and some 
constant C. This then yields Irn(A e . w (At)) > Ca £yUJ for x = 1 — £ 2 \n\ 2 /k 2 . Together with (52), 
we find 


dz 


C 


I {s 2 \k\ 2 /k^<l} 


dnlfuiz, k)\ < -3— ||/o(w,-)|| 22 (r2) . 

fc U," 


We turn now to the case z G (Lg, 0). We have 

1 


pI(z,k) = (al( L s, K)e ik “ x *’“W z / £2 + gg,(0, . 

Y ^£,Uj( K ) 


Owing (47), the first term of the r.h.s is direct and yields a control by C\\fo(oj, -)IIl 2 (r 2 )- For 
the second one, we write 


= m ^ _ a £ u( L s^) 

- rw w? j - : 

^£,ui(^') \j ^£,uj(k) 


(96) 


and obtain, thanks to (52), again a control by C\\fo(ui, •)||z. 2 (M 2 )- Consider finally the case 
z < Lg, for which 

0L(z,k) = e^“W L s- z V*fi,(Lg,K), z < Lg. 


The jump condition (42) yields p^(L s , At) = p^(Lg, At) — fo(u, At), which, together with (47) 
and (96), gives the expression 

pl{L~ s , At) = (i + e -2ik„\eAK)L S /£ 2 ^ + f u {Q^) e -i^sAK)L s /£\ 


Using again (52), we then proceed as in the case z > L and obtain the same estimate. Putting 
together all previous estimates ends the proof of the lemma. □ 
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Owing the previous lemma, we can proceed to the proof. Let v = P — p, which satisfies 


Av —“2 (l + £ s F(y, 7 ) 1 ( 0 ,L)( z )^dt v = ict £ p, (t, z, x ) € (0, + 00 ) x R x l 2 , 

equipped with u(0, •) = dtv( 0, •) = 0. Since V is uniformly bounded by a deterministic 
constant, and p £ C°((0, + 00 ), L 2 (M 3 )) according to Lemma 10.1, it is a classical problem to 
construct solutions to the above equation which satisfies the energy conservation relation 


1 d 

2 dt 


l|Vu(t)||i 2 (R 3 ) + 4 f 2 ( l + £ Sv ( 7 , j) 1 (o,L)^))\d t v(t,z,x)\ 2 dzdx] 

Cq «/M c c / 

= Re ia £ / pit, z, x)dtv(t, z, x)dxdz ) . 

V Jr 3 / 


After integration and the use of the Cauchy-Schwarz and Young inequalities, it follows that 

]|Vu(t)||i2 (R 3) + IldMOIli^RS) < Ca e J Q M s )Wl 3 (R 3 ) ds + C J Q ll^(s)||i a(R 3). 

Since v(t = 0, •) = 0, we can use the Poincare inequality (w.r.t. t), which together with the 
Gronwall lemma yield 

IKOIIhRrs) + 11^^)1112^3) < Cal I IIP( s )lli 2 (K3)ds- 


In order to apply Lemma 10.1 we notice that 

IIp(s)IIl 2 (K3) = £\\p(s,£-)\\ L 2 {R3) < £ \\p £ J\ L *(R 3 )duj. 

Above, we used the fact that p e u has compact support according to the aforementioned lemma 
since /o does. Standard Sobolev embeddings then yield 

sup ||'y(i, z, •)Hl 2 (ir 2 ) < Cal /2 , 

(t,z)e(0,T)xR 


which concludes the proof after rescaling x by ex. 
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